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Abstract 

An approach for synthesizing buckling results for, 
and behavior of, thin balanced and unbalanced symmetric 
laminates that are subjected to uniform heating or cooling 
and elastically restrained against thermal expansion or 
contraction is presented. This approach uses a nondimen- 
sional analysis for infinitely long, flexurally anisotropic 
plates that are subjected to combined mechanical loads 
and is based on useful nondimensional parameters. In ad- 
dition, stiffness-weighted laminate thermal-expansion pa- 
rameters and compliance coefficients are derived that are 
used to determine critical temperatures in terms of physi- 
cally intuitive mechanical-buckling coefficients. The ef- 
fects of membrane orthotropy and membrane anisotropy 
are included in the general formulation. Many results are 
presented for some common laminates that are intended to 
facilitate a structural designer’s transition to the use of ge- 
neric buckling design curves. Several curves that illus- 
trate the fundamental parameters used in the analysis are 
presented, for nine contemporary material systems, that 
provide physical insight into the buckling response in ad- 
dition to providing useful design data. Examples are pre- 
sented that demonstrate the use of generic design curves. 
The analysis approach and generic results indicate the ef- 
fects and characteristics of elastically restrained laminate 
thermal expansion or contraction, membrane orthotropy 
and anisotropy, and flexural orthotropy and anisotropy in 
a very general and unifying manner. 

Introduction 

Buckling behavior of laminated plates that are sub- 
jected to combined mechanical and thermal loads is an im- 
portant consideration in the design of advanced high- 
speed aerospace vehicles. The sizing of many structural 
subcomponents of these vehicles is often determined by 
stability constraints. One element that is of practical im- 
portance in structural design is the long rectangular plate. 
These plates commonly appear as elements of stiffened 
panels that are used for wing structures, and as semimono- 
coque shell segments that are used for fuselage and launch 
vehicle structures. Buckling results for infinitely long 
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plates are important because they often provide a practi- 
cal estimate of the behavior of finite-length rectangular 
plates, and they provide information that is useful in ex- 
plaining the behavior of these finite-length plates. More- 
over, knowledge of the behavior of infinitely long plates 
can provide insight into the buckling behavior of more 
complex structures such as stiffened panels. 

An important type ot long plate that appears as an 
element of advanced composite structures is the symmet- 
rically laminated plate. In the present study, the term, 
"symmetrically laminated," refers to plates in which ev- 
ery lamina, or ply, above the plate midplane has a corre- 
sponding lamina located at the same distance below the 
plate midplane, with the same thickness, material proper- 
ties, and fiber orientation. Symmetrically laminated 
plates are essentially flat after the manufacturing process 
and exhibit flat prebuckling deformation states, which is 
desirable for many applications. More importantly, the 
amenability of these plates to structural tailoring pro- 
vides symmetrically laminated plates with a significant 
potential for reducing the weight of aerospace vehicles or 
for meeting special performance requirements. Thus, 
understanding the mechanical and thermal buckling be- 
havior of symmetrically laminated plates is an important 
part of the search for ways to exploit plate orthotropy and 
anisotropy to reduce structural weight or to fulfill a spe- 
cial design requirement. 

In many practical cases, symmetrically laminated 
plates exhibit specially orthotropic behavior. However, 
in some cases, such as thin-walled [±45] s laminates that 
are candidates for spacecraft applications, these plates 
exhibit anisotropy in the form of material-induced cou- 
pling between pure bending and twisting deformations. 
This coupling is referred to herein as flexural anisotropy 
and it generally yields buckling modes that are skewed in 
appearance (see Fig. 1). Unbalanced, symmetrically 
laminated plates are also being investigated for special- 
purpose uses in aerospace structures. These laminated 
plates exhibit anisotropy in the form of material-induced 
coupling between pure inplane dilatation and inplane 
shear deformations in addition to flexural anisotropy. 
This coupling is referred to herein as membrane anisot- 
ropy and it generally yields combined inplane stress 
states for simple loadings like uniform edge compression 
when inplane displacement constraints are imposed on 
one or more edges of a plate. For example, when the edg- 
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es of an unbalanced, symmetrically laminated plate, such 
as a [+45 2 /0/90 ] s laminate, are totally restrained against 
thermal expansion and contraction, that is caused by uni- 
form heating or cooling, inplane shear stresses are devel- 
oped in addition to the usual compressive stresses that 
are often present in balanced laminates. These kinemati- 
cally induced shear stresses can be relatively large in 
magnitude, compared to the direct compressive stresses, 
and as a result can affect greatly the buckling behavior of 
the plate and yield buckling modes that are skewed in ap- 
pearance. 

The effects of flexural orthotropy and flexural 
anisotropy on the buckling behavior of long rectangular 
plates that are subjected to single and combined loading 
conditions are becoming better understood. For exam- 
ple, in-depth parametric studies that show the effects of 
flexural orthotropy and flexural anisotropy on the buck- 
ling behavior of long plates that are subjected to com- 
pression, shear, pure inplane bending, and various 
combinations of these loads have been presented in Refs. 

1 through 3. The results presented in these references in- 
dicate that the importance of flexural anisotropy on the 
buckling resistance of long plates varies with the magni- 
tude and type of the combined loading condition. Simi- 
lar results for plates loaded by uniform shear and a 
general linear distribution of axial load across the plate 
width have also been presented in Ref. 4. In a similar 
manner, the effects of membrane orthotropy and mem- 
brane anisotropy on the buckling behavior of long rect- 
angular plates that are restrained against axial thermal 
expansion or contraction and subjected to uniform heat- 
ing or cooling and mechanical loads have been presented 
in Refs. 5 and 6. Likewise, similar results for plates that 
are fully restrained against thermal expansion and con- 
traction and subjected to uniform heating or cooling have 
been presented in Ref. 7. This work has provided a bet- 
ter understanding of the load interaction effects of bal- 
anced and unbalanced, symmetrically laminated plates 
that are subjected to mechanical loads and restrained 
against thermal expansion or contraction. 

The effects of membrane orthotropy and anisotro- 
py, and the effects of flexural orthotropy and anisotropy, 
on the buckling behavior of long rectangular plates that 
are restrained against thermal expansion and contraction 
and subjected to uniform heating or cooling are becom- 
ing better understood. However, relatively little work 
has been done that addresses, in a broad way, the effects 
of a compliant, elastic restraining medium on the thermal 
buckling of plates that are restrained against expansion 
and contraction. This problem is important because, in 
reality, it is impossible to completely restrain a plate 
from thermal expansion and contraction, and in many in- 
stances the restraining medium undergoes thermal defor- 
mations. Thus, one objective of the present study is to 


present an analytical approach that complements the 
physically intuitive buckling analysis presented in Ret. 7 
for anisotropic plates that are fully restrained against 
thermal expansion or contraction and subjected to uni- 
form heating or cooling. Another objective is to indicate 
the effects of a compliant, elastic restraining medium in 
a very general manner. Towards these objectives, an in- 
tuitive buckling analysis for symmetrically laminated 
plates, that are elastically restrained against thermal ex- 
pansion and subjected to uniform heating or cooling, is 
presented that follows the analysis approach presented in 
Ref. 7, but includes the effects of the restraining medium 
in a very general manner. To achieve this objective, the 
buckling analysis is formulated directly in terms of buck- 
ling coefficients for the known, mechanically equivalent 
loads, stiffness-weighted laminate thermal expansion pa- 
rameters, and compliance coefficients for the restraining 
medium. Emphasis is placed on the properties ol the 
compliance coefficients for the restraining medium be- 
cause once known, they can be used with the generic re- 
sults presented in Ref. 7 for anisotropic plates that are 
fully restrained against thermal expansion and contrac- 
tion to determine the buckling behavior for a very broad 
range of laminate constructions and restraining media. 
Thus, another objective of the proposed paper is to 
present values of the compliance coefficients for the re- 
straining medium for a wide range of laminate stacking 
sequences and material systems. 

To accomplish the objectives of the present study, 
the analysis that is used to formulate thermal buckling 
problem in a generic nondimensionat form is described 
first. This description includes the nondimensionat 
buckling analysis for infinitely long symmetrically lam- 
inated plates that are subjected to combined mechanical 
loads, and the formulation for the thermally induced me- 
chanical loads and critical temperature change. Then, 
examples of compliance coefficients for a laboratory- 
type elastic restraining medium are presented. Next, spe- 
cific results are presented for several common laminates 
with the two long edges clamped or simply supported 
and all edges elastically restrained against inplane move- 
ment. Generic results are then presented that illustrate 
the fundamental parameters that are used in the analysis 
for nine contemporary material systems. Lastly, exam- 
ples are presented that use the generic buckling-design 
curves presented herein and in Ref. 7, that are applicable 
to a very wide range of laminate constructions, to dem- 
onstrate the analytical procedure and illustrate the effects 
of a restraining medium on the buckling behavior. 

Analysis Description 

In preparing generic design charts for buckling of a 
single flat thin plate, a special-purpose analysis is often 


2 

American Institute of Aeronautics and Astronautics 



preferred over a general-purpose analysis code, such as a 
finite-element code, because of the cost and effort that is 
usually involved in generating a large number of results 
with a general-purpose code. The results presented in the 
present study were obtained by using such a special-pur- 
pose buckling analysis that is based on classical laminat- 
ed-plate theory. The analysis details are lengthy; hence, 
only a brief description of the buckling analysis is pre- 
sented herein. First, the buckling analysis for long plates 
that are subjected to a general set of mechanical loads is 
described. Then, the mechanical loads that are induced 
by elastically restrained thermal expansion and contrac- 
tion and that are used in the buckling analysis are de- 
rived, and an expression for the critical temperature 
change is presented that is in terms of the corresponding 
crtitical loading parameter and mechanical-buckling co- 
efficients. 


Symmetrically laminated plates can have many dif- 
ferent constructions because of the wide variety of mate- 
rial systems, fiber orientations, and stacking sequences 
that can be selected to construct a laminate. A way of 
coping with the large number of choices for laminate 
constructions is to use convenient nondimensional pa- 
rameters in order to understand overall behavioral trends 
and sensitivities of the structural behavior to changes in 
laminate construction. The buckling analysis used in the 
present paper is based on classical laminated-plate theo- 
ry and the classical Rayleigh-Ritz method, and is derived 
explicitly in terms of the nondimensional parameters de- 
fined in Refs. 1-8. This approach was motivated by the 
need for generic (independent of a specific laminate con- 
struction) parametric results for composite-plate buck- 
ling behavior that are expressed in terms of the minimum 
number of independent parameters needed to fully char- 
acterize the behavior, and that indicate the overall trends 
and sensitivity of the results to changes in the parame- 
ters. The nondimensional parameters that were used to 
formulate the buckling analysis are given by 
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where b is the plate width and X is the half-wave 
length of the buckle pattern of an infinitely long plate 
(see Fig. 1 ). The subscripted D-terms are the bending 


stiffnesses of classical laminated-plate theory. The 
parameters a„, and (3 characterize the flexural orthot- 
ropy, and the parameters y and 5 characterize the flex- 
ural anisotropy. 

The mechanical loading conditions that are includ- 
ed in the buckling analysis are uniform transverse ten- 
sion or compression, uniform shear, and a general linear 
distribution of axial load across the plate width, as de- 
picted in Fig. 1. Typically, an axial stress resultant dis- 
tribution is partitioned into a uniform part and a pure 
bending part. However, this representation is not unique. 
The longitudinal stress resultant N x is partitioned in the 
analysis into a uniform tension or compression part and 
a linearly varying part that corresponds to eccentric in- 
plane bending loads. This partitioning is given by 

N l = N«-N b [e 0 + (e,-e 0 )Ti] (5) 

where N xc denotes the intensity of the constant-valued 
tension or compression part of the load, and the term 
containing N b defines the intensity of the eccentric 
inplane bending load distribution. The symbols £<, and 
£j define the distribution of the inplane bending load, 
and the symbol r| is the nondimensional coordinate 
given by r| = y/b. This particular way of partitioning 
the longitudinal stress resultant was used for conve- 
nience by eliminating the need to calculate the uniform 
and pure bending parts of an axial stress resultant distri- 
bution prior to performing a buckling analysis. 

The analysis is based on a general formulation that 
includes combined destabilizing loads that are propor- 
tional to a positive-valued loading parameter p that is 
increased until buckling occurs, and independent suberit- 
ical combined loads that remain fixed at a specified load 
level below the value of the buckling load. Herein, the 
term "subcritical load” is defined as any load that does 
not cause buckling to occur. In practice, the subcritical 
loads are applied to a plate prior to, and independent of, 
the destabilizing loads with an intensity below that which 
will cause the plate to buckle. Then, with the subcritical 
loads fixed, the active, destabilizing loads are applied by 
increasing the magnitude of the loading parameter until 
buckling occurs. This approach permits certain types of 
combined-load interaction to be investigated in a direct 
and convenient manner. For example, in analyzing the 
stability of an aircraft fuselage, the nondestabilizing 
transverse tension load in a fuselage panel that is caused 
by cabin pressurization can be considered to remain con- 
stant and, as a result, it can be represented as a passive, 
subcritical load. The combined shear, compression, and 
inplane bending loads that are caused by flight maneu- 
vers can vary and cause buckling and, as a result, they 
can be represented as active, destabilizing loads. 

The distinction between the active, destabilizing 
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and passive subcritical loading systems is implemented 
in the buckling analysis by partitioning the prebuckling 
stress resultants as follows 


N xe =-N' xl +N * l2 (6) 

N y = — N y , + N y2 (7) 

N xy = N xyl + N xy2 (8) 

N b = N bl +N b2 (9) 


where the stress resultants with the subscript 1 are the 
destabilizing loads, and those with the subscript 2 are 
the subcritical loads. The sign convention used herein 
for positive values of these stress resultants is shown in 
Fig. 1. In particular, positive values of the general linear 
edge stress distribution parameters N hl , N b2 , e {) , and 
correspond to compression loads. Negative values of 
N hl and N b2 , or negative values of either or e,, yield 
linearly varying stress distributions that include tension. 
Depictions of a variety of inplane bending load distribu- 
tions are given in Ref. 4. The two normal stress result- 
ants of the system of destabilizing loads, N xl and N y] , 
are defined to be positive-valued for compression loads. 
This convention results in positive eigenvalues being 
used to indicate instability due to uniform compression 
loads. 

The buckling analysis includes several nondimen- 
sional stress resultants associated with Eqs. (6) through 
(9). These dimensionless stress resultants are given by 
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the specific form (relative contributions of the load com- 
ponents) of a given system of destabilizing loads. Typi- 
cally, the dominant load factor is assigned a value of 1 
and all others are given as positive or negative fractions. 

Nondimensional buckling coefficients that are used 
herein are given by the values of the dimensionless stress 
resultants of the system of destabilizing loads at the onset 
of buckling; i.e., 
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where the quantities enclosed in the parentheses with 
the subscript “cr M are critical values that correspond to 

buckling and p cr is the magnitude of the loading param- 
eter at buckling. Positive values of the coefficients K x 
and Ky correspond to uniform compression loads, and 
the coefficient K s corresponds to uniform positive shear. 
The direction of a positive shear stress resultant that acts 
on a plate is shown in Fig. 1. The coefficient K b corre- 
sponds to the specific inplane bending load distribution 
defined by the selected values of the parameters e (> and 
£, (see Fig. I). 

The mathematical expression used in the variation- 
al analysis to represent the general off-center and skewed 
buckle pattern is given by 
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where the subscript j takes on the values of 1 and 2. In 
addition, the destabilizing loads are expressed in terms 

of the loading parameter p in the analysis by 
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where L] through L 4 are load factors that determine 


w N (£,t|)= £ (A m sin^ + B m cosj^)<l> m (ri) ( 22 ) 

where £ = \/X and r| = y/b are nondimensional coor- 
dinates, w N is the out-of-plane displacement field, and 

A m and B m are the unknown displacement ampli- 
tudes. In accordance with the Rayleigh-Ritz method, 

the basis functions O m (T|) are required to satisfy the 
kinematic boundary conditions on the plate edges at r\ = 
0 and 1 . For the simply supported plates, the basis 
functions used in the analysis are given by 

O m (T|) - sin mTcrj (23) 

for values of m = 1, 2, 3, ..., N. Similarly, for the 
clamped plates, the basis functions are given by 
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form of the membrane constitutive equations for sym- 
<J> m (r|) = cos(m- 1 )7tn - cos(m+ 1 ) 7 tri (24) metrically laminated plates is given by 


For both boundary conditions, the two long edges of a 
plate are free to move in-plane. 

Algebraic equations that govern the buckling be- 
havior of infinitely long plates are obtained by substitut- 
ing the series expansion for the buckling mode given by 
Eq. (22) into a nondimensionalized form of the second 
variation of the total potential energy and then comput- 
ing the integrals appearing in the nondimensional second 
variation in closed form. The resulting equations consti- 
tute a generalized eigenvalue problem that depends on 

the aspect ratio of the buckle pattern )Jb (see Fig. 1 ) and 
the nondimensional parameters and nondimensional 
stress resultants defined herein. The smallest eigenvalue 
of the problem corresponds to buckling and is found by 

specifying a value of X/b and solving the corresponding 
generalized eigenvalue problem for its smallest eigen- 
value. This process is repeated for successive values of 

X/b until the overall smallest eigenvalue is found. 

Results that were obtained from the analysis de- 
scribed herein for uniform compression, uniform shear, 
pure inplane bending (given by £<,= -1 and E x = 1), and 
various combinations of these mechanical loads have 
been compared with other results for isotropic, orthotro- 
pic, and anisotropic plates that were obtained by using 
other analysis methods. These comparisons are dis- 
cussed in Refs. 1-3, and in every case the results de- 
scribed herein were found to be in good agreement with 
those obtained from other analyses. Likewise, results 
were obtained for isotropic and specially orthotropic 
plates that are subjected to a general linear distribution of 
axial load across the plate width and compared with re- 
sults that were obtained by seven different authors (see 
Ref. 4). In every case, the agreement was good. 

Prebuckling Stress Resultants 

Symmetrically laminated plates that are subjected 
to uniform temperature fields and restrained against ther- 
mal expansion and contraction may develop internal me- 
chanical loads that can cause buckling. These induced 
mechanical loads enter the analysis through the mem- 
brane constitutive equations. The standard form of these 
membrane constitutive equations for thin plates, that is 
based on classical laminated-plate theory, is found in 
Refs. 9 and 10 and is expressed in terms of membrane 
stiffness coefficients and fictitious thermal stress result- 
ants. Following Ref. 7, the membrane constitutive equa- 
tions are used in the present study in an inverted form 
that uses the overall laminate coefficients of thermal ex- 
pansion and the membrane compliance coefficients. This 



a n a i2 a i6 
a i2 a 22 a 26 
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where u(x,y) and v(x,y) are the prebuckling, inplane dis- 
placements in the x- and y-coordinate directions (see 

Fig. 2), respectively; a x , a y , and a xy are the overall 
laminate coefficients of thermal expansion; the sub- 
scripted a-terms are the plate membrane compliance 

coefficients; 0 O is the magnitude of the uniform tem- 
perature change from a predetermined stress- and strain- 
free reference state; and commas followed by a sub- 
script denote partial differentiation with respect to the 
coordinate associated with the subscript. For restrained 
thermal expansion and contraction problems, the plates 
are assumed to be supported and loaded such that the 
prebuckling stress field is uniform, that is, homoge- 
neous. With this assumption, a compatible displacement 
field is obtained directly by integrating Eqs. (25). This 
integration yields 

u(x,y) = |a„ N x + a l2 N y + a, 6 N xy + oq ©„jx + g,y + g 2 (26) 

v(x,y) = |a, 2 N„ + a 22 N y + a 26 N xy + a, 0„jy + g,x + g 4 (27) 

and 

gi + g} = a i6 N x + a 26 N y + a 66 N xy + a xy 0 (1 (28) 

Equations (25)-(28) can be used to determine the 
thermally induced mechanical loadings for several prob- 
lems of practical interest. The problems consist of plates 
restrained against axial thermal expansion and contrac- 
tion (see Ref. 5 or 6), plates restrained against transverse 
thermal expansion or contraction (y-coordinate direc- 
tion), and plates completely restrained against thermal 
expansion and contraction (see Ref. 7). In the present 
study, however, symmetrically laminated plates that are 
elastically restrained against thermal expansion and con- 
traction are considered. 

All the stress resultants in a given plate are in- 
duced by the elastically restrained thermal expansion or 
contraction and, when considered together, may be de- 
stabilizing. The induced stress resultants are proportion- 
al to the strains caused by expansions or contractions and 
shearing deformations of the plate and the resistance pro- 
vided by the restraining medium. Typically, the elastic 
resistance of a restraining medium is simulated with lin- 
ear springs and expressed in terms of the corresponding 
spring stiffnesses. In the present study, the elastic resis- 
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tance of a restraining medium is represented approxi- 
mately, and in a general manner, in terms of the overall 
compliance and thermal expansion coefficients of a ho- 
mogeneous, anisotropic restraining medium. For this 
type of representation, the prebuckling stress state is as- 
sumed to be uniform, for all practical purposes. When 
the compliance of the restraining medium is relatively 
small (close to a rigid medium), the induced prebuckling 
stress state is expected to be nearly uniform and the ap- 
proach presented herein for calculating buckling loads is 
expected to be a reasonable approximation. For more 
compliant restraining media, nonuniformities in the pre- 
buckling stress state are expected to increase, and thus 
the approach presented herein for calculating buckling 
loads must be used with caution and considerable engi- 
neering judgement. However, it is apparent that there are 
several types of elastic support conditions that can be 
modeled accurately with the approach presented herein. 
One type of elastic support condition, that may be useful 
in laboratory structural testing, is described herein after 
the general formulation of the induced prebuckling stress 
state and buckling analysis are presented. 

In the genera] formulation presented herein, the 
thermally induced inplane stress resultants are deter- 
mined by effectively enforcing displacement compatibil- 
ity and Newton’s law of action and reaction at the 
interfaces between the plate and the restraining medium. 
In particular, the stresses in the plate and the restraining 
medium are assumed to be, for the most part, uniform 
and to have the same constant values. Moreover, elonga- 
tion and shear of the plate are assumed to be balanced by 
contraction and shear of the restraining medium such that 
the total, average strain in the plate is the negative of the 
total, average strain in the part of the restraining medium 
providing the resistance. For a restraining medium like 
the one shown in Fig. 2, but of effectively infinite planar 
extent, the stresses produced by expansion or contraction 
and shear are self equilibrated and, as a result, the elastic 
resistance provided can be idealized, or approximated, as 
the result of a uniform stress and deformation state acting 
over a rectangular region surrounding the plate and with 
finite dimensions. The finite dimensions could be deter- 
mined by appyling St. Venant’s principle to determine 
the decay lengths of the self-equilibrated stresses. 

Following the reasoning of the previous supposi- 
tions, the positive-valued (total) strains in the restraining 
medium are obtained from the equivalent of Eq. (25) for 
the restaining medium. Then, the strains in the restrain- 
ing medium are replaced in the equivalent of Eq. (25) 
with the negative of the corresponding strains in the 
plate, given by the left-hand side of Eq. (25), to fulfill 
displacement compatibilty in an overall manner. With 
this approach, the positive-valued (total) strains in the 
plate, given by the left-hand side of Eq. (25), are ex- 


pressed in terms of the properties of the restraining me- 
dium as 
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where the superscript R denotes a quantity associated 
with the restraining medium. It is reiterated, that the 
negative signs in Eq. (29) come from the supposition 
that a positive, expansional strain for the plate in the x- 
direction produces compressive stress in the x-direction, 
that corresponds to a negative value for N x . In addition, 
the Poisson effect produces a tensile stress in the y- 
direction, that corresponds to a positive value for N y , A 
similar response is noted for a positive, expansional 
strain for the plate in the y-direction. For a positive 
shearing deformation, a negative shearing stress result- 
ant acts on the restraining medium. 

Enforcing compatibility at the interface between 
the plate and the restraining medium is manifested by 
eliminating the strains from Eqs. (25) and (29), which 
gives 


(a„+a" ) (a tl +af 3 ) (a 16 +ai) 
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Next, because all the subcritical loads are zero-valued in 
the absence of applied mechanical loads, N x (x,y) = 

-Nxi , N y (x,y) = -N yl , and N xy (x,y) = N xyl (see Fig. 2). 
This substitution yields 

(a lk +ar.) (a 12 +a*) (a l6 +a? 6 ) 

(a, 2 +a* 2 ) (a 22 +a 21 ) (a 26 +a 2h ) 

(a 16 +af 6 ) (a 2A +a*) (a w +a^) 

The thermally induced stress resultants for this problem 
are obtained by solving matrix Eq. (31) for Nil , N yl , and 
N The solution is given symbolically by 

N:, = [£ n K + a?) + A l2 (a y + <) + A „K + <) ]0„ (32) 

N„ = [ A 12 (ot K + c£) + A 22 (a y + ct?) + A + <) ] 0 () (33) 

N*=-[X w (a. + o£) +A 26 (a y + a y R ) +A«(a B , + <)]0« (34) 


N', \ 

I 0,+aJ 

Nyl } = < 

a y +a; 

-nJ 

\ a,y+< 


K (31) 


where the subscripted A -terms are elements of the 
inverse of the matrix in the left-hand side of matrix Eq. 
(31). When the restraining medium is totally rigid, the 
compliances of the restraining medium in Eq. (31) are 
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zero valued and the subscripted A -terms reduce to the 
usual membrane stiffnesses of laminated plate theory, 
which are denoted herein by subscripted A-terms. 
Because the elements of the matrix in the left-hand side 
of Eq. (31) typically have greater magnitude than the 
corresponding elements for only the plate, and because 
the elements of the inverse of a matrix are like their 
corresponding reciprocals, it follows that the sub- 
scripted A -terms typically have smaller magnitude than 
the corresponding subscripted A-terms. 

Critical Temperature Change 

Equations (32)-(34) define a combined loading 
state that is induced by elastically restrained thermal ex- 
pansion and contraction. More specifically, the equa- 
tions indicate that the induced loads depend on the 
compliance of the restraining medium and the relative 
overall thermal expansion or contraction of the plate and 
the restraining medium. In addition, these equations 
show that each of the thermally induced mechanical 
loads generally depends on all three laminate coeffi- 
cients of thermal expansion, and that positive values for 
the stress resultants are possible even for negative values 
of © () (uniform cooling) and vice versa. For example, a 
laminate could have a negative value of a x and still have 

a positive value of N*i (axial compression). Thus, the 
signs of 0 O and the bracketed quantities in Eqs. (32)-(34) 
must be considered in formulating the buckling problem. 

The buckling problem is posed by first substituting 
Eqs. (32)-(34) into Eqs. (10)-(12), respectively, to obtain 
expressions for the nondimensional stress resultants that 
are used to characterize the thermally induced mechani- 
cal loads. Following Ref. 7, the nondimensional stress 
resultants are expressed in terms of stiffness-weighted 

laminate thermal-expansion parameters denoted by &i, 

<5t 2 , and <5t 3 , and in terms of nondimensional compliance 
coefficients denoted by C p C 2 , and C 3 . These expres- 
sions, with the use of Eqs. (14)-(16), are given by 


where 


„ c 1 2b pk __ i k 

7tV (1 +c,)° 0 ,p 

(35) 

n 2 t 2 (l +C 2 )°° 2P 

(36) 

1 2b pi T pc 

n *’ , “ n 2 t 2 (1 + C,)° 0_L5p 

(37) 

_ __ t (A i | 0 t x + A l2 ct y + A, ft ot xy ) 
' 12VD..D* 

(38) 


C,= 


c 2 = 




_ r(A 12 a„ + A 22 a y + A 26 a, y ) 
® 2 ~ 12D,, 

t 2 (A l 6 a, + A M a y + A^a,,) 

3 12 (D„ Dj 2 ) 1m 

A„a, + A i 2 a, + A l 6 a t , 

A lt (a, + aj) + A 12 (a y + a y ) + A „(«» + <) 

A l 2 a,+ A^-t A Jft a, y 
A 12 (a, + a?) + A 22 (a y + a") + A 26 (a, y + O 

A l<s a, + A 2 <s a y + Ajx xy 

A w ( a, + a?) + A 2 fc (a y + a y R ) + A<*(a, y + a* ) 


(39) 

(40) 

(41) 

(42) 

(43) 


It is important to point out that because the sub- 
scripted A-terms are typically less than the correspond- 
ing subscripted A-terms, that typically Cj > 0, C 2 > 0, and 
C 3 >0 when the thermal expansion or contraction of the 
restraining medium is negligible. Moreover, for a totally 
rigid restraining medium with negligible thermal expan- 
sion or contraction, C, = C 2 = C 3 = 0. However, in gen- 
eral, the compliance factors C,, C 2 , and C 3 may be 
negative valued when the restraining medium has a high- 
er relative thermal expansion or contraction than the 
plate. Thus, the relative stiffness and relative thermal ex- 
pansion and contraction can significantly change the na- 
ture of the thermally induced mechanical loads, as 
expected. 

Equations ( 1 0)-( 1 2) and (35)-(37) indicate that 
Nxi, N |, andN xyl are positive-valued (see Fig. 2b) when 

&,/(! + Ci), ft 2 /(l + C 2 ), and a 3 /(l + C 3 ) are positive- 
valued, respectively, and when 0 O is positive-valued. 
Similarly, Nh, N yl , and N xyl are negative-valued when 

&i/(l + Ci), Ol 2 /{\ + C 2 )» and & 3 /(l + C 3 ) are positive- 
valued, respectively, and 0 () is negative-valued, or when 

<5ti/(l +C,), ft 2 /(l +C 2 ), and a 3 /(l + C 3 ) are negative- 
valued, respectively, and 0 O is positive-valued. Next, 
Eqs. (35)-(37) are substituted into Eqs. (18)-(20) to ob- 
tain the relationships between the mechanical-buckling 
coefficients, load factors, and the critical temperature 

change 0o; that is. 


K --("“,- = W(T^ e ‘ = L||> " 

= ( n yl ) tT = ^ ^ + c 2 ) 0 <> = L > P« 


(44) 

(45) 

(46) 
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where the critical eigenvalue p cr = p LT (P, Y» 5, L 1? L 2 , L 3 ) 
for a given set of flexural (bending) boundary conditions 
(e.g., simply supported and clamped edges). Eqs. (44)- 
(46) reduce to the corresponding equations presented in 
Ref. 7 for plates that are fully restrained against thermal 
expansion and contraction (C, = C 2 = C 3 = 0). For an 

isotropic material, - 0 and ft, = ft 2 = <x( 1 + v) , where 
a is the coefficient of thermal expansion and v is Pois- 
son’s ratio of a homogeneous isotropic material. 

The next step in posing the buckling problem is to 
define the load factors L,, L 2 , and L 3 that appear in Eqs. 

(35) -(37) and (44)-(46). It is important to reiterate that 

positive, negative, and zero values for fti0<>/(l +Ci) 
correspond to positive, negative, and zero values for Nxi» 
respectively (see Fig. 2b). Similarly, positive, negative, 

and zero values for ft 2 0o/(l + C 2 ) correspond to posi- 
tive, negative, and zero values for N yl , respectively; and 

positive, negative, and zero values for ft 3 0o/(l + C 3 ) 
correspond to positive, negative, and zero values for 
N p respectively. To define the load factors properly, 

the signs of N£i, N yl , and N xyl must be considered. Spe- 
cifically, the load factors must be defined such that pos- 
itive values of L,, L 2 , and L 3 correspond to positive 

values of Nlu N yl , and N xyl , respectively. Moreover, 
both positive (heating) and negative (cooling) values of 
0 O must be considered. These requirements lead to six 
cases that must be considered in formulating the buck- 
ling analysis; that is, the cases for which 

fii0„/(l +Ci) > 0, ft!0o/( 1+Ci)<0, ft, =0 with 
ft 2 0«/(l+C 2 )> 0, ft, = 0 with ft 2 0 o/(l+C 2 )< 0, 
ftj = ft 2 = 0 with ft 3 0o/( 1 + C 3 ) > 0, and ft, = ft 2 = 0 

with ft 3 0 ( /( 1 + C 3 ) < 0. The values of the mechanical- 
load factors L,, L 2 , and L 3 for each of these cases are pre- 
sented subsequently. 

Load factors for case 1 . For the case when 
ft,0o/(l + Ci) > 0 , Nxi > 0 and Lj = 1 is appropriate 
(axial compression). The values for the other two load 
factors that are needed to completely define the preb- 
uckling stress state are obtained by dividing Eqs. (19) 
and (20) by Eq. (18), with L, = 1, or by dividing Eqs. 

(36) and (37) by Eq. (35). This step yields 



r 1 +C,\ A, 6 cc, + A 26 a y + A M ff g |D„\»" 
J + cj A u a % + A 12 cty + A D 2 J 


For an isotropic plate, these expressions reduce to 
L 2 = (1 + C 1 )/( 1 + C 2 ) and L 3 = 0. With L, = 1 and 
L 2 and L 3 defined by Eqs. (47) and (48), the relation- 
ship between the critical value of the mechanical load- 
ing parameter p cr and the critical temperature change 
0o is determined by Eq. (44); that is. 


K - 12b 2 a, 

* 7tV (1 + c,) 


= Pcr 


(49) 


where p cr = p cr (P, Y> 8, L ,, L 2 , L 3 ) for a given set of 
flexural boundary conditions. It is important to point 
out that Eq. (49) yields positive and negative values for 

0o for positive and negative values of ft,/(l +C,), 
respectively. Moreover, it is important to reiterate that 

the relationship between p cr and the corresponding 
mechanical buckling coefficients K x , K y , and K s is given 
by Eqs. (49), (45), and (46), respectively. 

Load factors for case 2 . For the case when 
ft i0o/( 1 + Ci) < 0, N*i < 0 and L, = -1 is appropriate 
(axial tension). Like for the previous case, the values for 
the other two load factors that are needed to completely 
define the prebuckling stress state are obtained by divid- 
ing Eqs. (19) and (20) by Eq. (18), but with L, = -1 , or by 
dividing Eqs. (36) and (37) by Eq. (35). This step yields 



1 +C, ) A„.a, + A ;t a, + A M a„ /D„\ w 
I +cjA„a I + A, i a y + A, 6 a„(Dj 



1 +C,\ A|,a, + A„a, + A ; .a, y / D l | \ 1 ' 1 

1 + C J A , ,a. + A l2 a, + A ,„a„ ( D„ j 


(47) 


With L, = -1 and L 2 and L 3 defined by Eqs. (50) and 
(51), the relationship between the critical value of the 
mechanical loading parameter p cr and the critical tem- 
perature change 0« is again determined by Eq. (44); 
that is, 
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K = 


12b 2 a, _ 
TtV (i + c ,) 0 


-Per 


(52a) 


In contrast to the previous case, Eq. (52a) yields positive 

values for ©o for negative values of fti/(l +Cj), and 
vice versa. For laminates with L 3 = 0 (balanced lami- 
nates) and L 2 < 0 (no transverse compression), no desta- 
bilizing compression or shear loads are present and 
buckling cannot occur because the plate is in a state of 
either axial or biaxial tension. In contrast, when L 3 = 0 
and L 2 > 0, a plate is subjected to axial tension and 
transverse compression. Figures 29-31 of Ref. 1 indi- 
cate that an infinitely long plate buckles as a wide col- 
umn for this type of loading, and that the buckling 
coefficient Ky = 1 and 4 for simply supported and 
clamped plates, respectively. With Ky known, Eqs. (45) 
and (52a) give 


K = 


12b 2 &i qct _ __ 

TcVfl+C,) °“ L 2 


(52b) 


It is important to point out that Eq. (54a) yields positive 
and negative values for ©o for positive and negative 

values of & 2 /(l + C 2 ), respectively. For laminates with 
L 3 = 0 (balanced laminates) a plate is subjected to only 
transverse compression. Thus, an infinitely long plate 
buckles as a wide column for this type of loading, and 
the buckling coefficient K y = 1 and 4 for simply sup- 
ported and clamped plates, respectively. With K y 
known, Eq. (45) gives 


12 b 2 <* 2 

n 2 t 2 (1 + C 2 ) 


©n = K v 


(54b) 


Load factors for case 4 . For the case when ft, = 0 

and & 2 0o/(l + C 2 ) < 0, N yl < 0 and L, = 0 and L 2 = 
-1 are appropriate (transverse tension). Like for the pre- 
vious case, the value for the load factor L 3 that is needed 
to completely define the prebuckling stress state is ob- 
tained by dividing Eq. (20) by Eq. (19), with L 2 = -1, or 
by dividing Eqs. (37) by Eq. (36). This step yields 


Load factors for case 3 . For a general symmetric 
laminate, the possibility exists that <5ti = 0 , which im- 
plies that Nxi =0. For this case, = 0 is appropriate 

and the sign of ft 2 0o/( 1 + C 2 ) must be considered in 
defining the nonzero load factors. In particular, for 
a 2 0 ( >/( 1 + C 2 ) ^ 0 , Ny| ^ 0 and Lj = 0 and L 2 — 1 are 
appropriate (transverse compression). The value for the 
load factor L 3 that is needed to completely define the 
prebuckling stress state is obtained by dividing Eq. (20) 
by Eq. (19), with L 2 = 1, or by dividing Eqs. (37) by Eq. 
(36). This step yields 



1 + C 2 A „a t + A 2b a y + A M a xy / p 22 •» 
1 + cj A^a, + A 22 a y + A 2 ft a iy \ D n 


The relationship between L 3 and the mechanical-buck- 
ling coefficients Ky and K s is shown in Figs. 24-27 of 
Ref. 1. With L, = 0, L 2 =1, and L 3 defined by Eq. (53), 
the relationship between the critical value of the 
mechanical loading parameter per and the critical tem- 
perature change 0o r is determined by Eq. (45); that is, 


K = 


12 b 2 6^2 pyCT 

nh 2 (1 +C 2 ) 0 


= P« 


(54a) 


T — _ 

- M 


N.„ /D 

"V- 

1 1 +c,ja, 

N* (d 

J " 

U+cJa, 


[1+C , 1 

l A „,a, + A K a, + A j D„ 


i 1 + C)j 

|A, ; a, + A^a, + A !( ,a >y \D„ 


Like the previous case, the relationship between L 3 and 
the mechanical-buckling coefficients Ky, and K s is shown 
in Figs. 24-27 of Ref. 1. With L, = 0, L 2 = -1, and L 3 
defined by Eq. (55), the relationship between the critical 
value of the mechanical loading parameter p cr and the 
critical temperature change 0? is determined by Eq. 
(45); that is, 


K = 


12 b 2 

n\ 2 (1 + C 2 ) 


0 " = -p cl 


(56) 


For this case, Eq. (56) yields positive values for 0o for 

negative values of <x 2 /(l + C 2 ), and vice versa. For lam- 
inates with L 3 = 0 (balanced laminates), no destabilizing 
compression or shear loads are present and buckling 
cannot occur because the plate is in a state of uniaxial, 
transverse tension. 

Load factors for case 5 . For the case with 
ft, =& 2 = 0 and & 3 0o/O + C 3 ) > 0 , N c *i = N yl =0 and 
N xy , > 0, which implies that Lj = L 2 = 0, and that L 3 = 1 
is appropriate (positive shear loading). With L, = 0, L 2 
= 0, and L 3 = 1, the relationship between the critical val- 
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ue of the mechanical loading parameter p cr and the crit- 
ical temperature change 0o r is determined by Eq. (46); 
that is, 


K s = 


12b 2 

nh 2 (1 +C 3 ) 


00= Pa 


(57) 


Again, it is important to point out that Eq. (57) yields 
positive and negative values for 0[> r for positive and 

negative values of <5t 3 /(l + C 3 ), respectively. In addi- 
tion, values of the mechanical-buckling coefficient K s 
for several laminates are given in Ref. 1 . 

Load factors for case 6 . The final case to consider 

is when 6ti = &2 = 0 and & 30 «/(l + C 3 ) < 0. For this 

case, N*i = N y i = 0 and N xy] < 0, which implies that 
= L 2 = 0, and that L 3 = - 1 is appropriate (negative shear 
loading). With L, = 0, L 2 = 0, and L 3 = -1, the relation- 
ship between the critical value of the mechanical loading 
parameter p cr and the critical temperature change ©o is 
again determined by Eq. (46); that is, 


K = 


12b 2 ^ 

K 2 1 2 (1 + Q) 


00= -Pa 


(58) 


In contrast to the previous case, Eq. (58) yields positive 
values for 0» for negative values of & 3 /(l + C 3 ) , and 
vice versa. 

Finally, it is important to mention that the approach 
used herein to define the prebuckling stress state and the 

critical temperature change 0o also applies for a more 
sophisticated plate theory, like a first-order transverse- 

shear deformation theory. For this theory, p cr would de- 
pend also upon additional nondimensional parameters 
that characterize the transverse-shear flexibility. Thus, 
the only difference in the results for the two plate bend- 
ing theories is the actual value of p cr that is used in Eqs. 
(44)-(46), for a given problem. It is also important to 
point out that p tr for an infinitely long plate does not de- 
pend on the parameter a^. This fact has been shown in 
Refs. 1-4. 


An idealized structural arrangement is shown in 
Fig. 3 in which a central balanced, symmetrically lami- 
nated plate is restrained against inplane movement by 
four generally different, balanced, symmetrically lami- 
nated plates. For this arrangement, the thermally in- 
duced loads consist of tension or compression loads, 
without shear. The plates that provide the elastic re- 
straint are number sequentially from one to four in the 
figure. Each of the four restraining plates is subjected to 
a generally different temperature change, denoted by 

0o° (k = 1, 2, 3, and 4), and is enclosed by an insulated 
boundary. In addition, the central plate is subjected to a 
temperature change, denoted by 0 O , where the super- 
script (5) has been dropped for simplicity. The restrain- 
ing plates are supported to allow tangential movement at 
the rigid frame shown in the figure and at the interfaces 
with the central plate. 

The compliance of the restraining medium is deter- 
mined by first defining the elongations of each of the five 
plates. In particular, the axial elongation of plates 1, 2, 

and 5 are denoted by 5 ( x l) , 8 X 2) , and 6 * 5) , respectively. 
Similarly, the transverse elongation of plates 3, 4, and 5 

are denoted by 8 y 3) , 8 y 4) , and 8 y respectively. To satisfy 
compatibility, the elongations must satisfy the following 
equations 

5 X U + 5 X 2) + S ( x 5) = 0 (59) 

5y 3) + 8y 4) + 5y ?) = 0 (60) 

because the surrounding frame is rigid. These compati- 
bility conditions state that the axial and transverse elon- 
gation or expansion of the central plate is balanced by 
the corresponding contraction of the restraining plates. 
These conditions are expressed in terms of strains by 

^e< n + a '-e< 2) + ei” = 0 (61) 

1 O) l(4) 

^- e o) + ^_4 4) + e^ 5) = 0 (62) 

where e ( x k) and e ( y k) (k = 1 , 2, 3, 4, and 5) are given by the 
corresponding elongation divided by the corresponding 
length or width. 

Next, the strains in each plate are expressed in 
terms of the stress resultants as follows 


Example of Compliance Coefficients 

Before performing the buckling analysis described 
herein, expressions for the compliance coefficients de- 
fined by Eqs. (41)-(43) are needed. To demonstrate the 
formulation of the present study, the example of a prac- 
tical, laboratory-type apparatus, that is depicted in Fig. 3, 
is presented in this section. 


e?» = a')VN'. k) + aiVNy 1 + ai kl 0» kl (63) 

e< k > = a)VN ,k> + a©N; k) + a< k) 0, ( !° (64) 

where the superscript (k) identifies the quantity for a 
given plate. These expressions are simplified by noting 

that N x 3) = N ( x 4) = 0 and N y n = N y 2) = 0 , and by applying 
Newton’s law of action and reaction between adjoining 
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plates, which gives Ni l) = = Nf } = N x and 

Ny 3) = Ny 4) = Ny } = N y . Substituting these expression 
for the stress resultants into Eqs. (63) and (64) and then 
the resulting expressions for the strains into Eqs. (61) 
and (62) gives 

(a n + + ^pafi’jNx + a, 2 N, = 

- (a x 0» + x a > l)0| > n + ^f-a?©®) (65) 

ai 2 N x + (a 22 + + ^-a$jN y = 

- (a y 0„ + ^a' 3 ©!, 3 ’ + ^a' 4 ©!?’) (66) 


compliance and thermal expansion or contraction of the 
central plate and the restraining plates. For Eq. (72), 
Eqs. (32) and (33) reduce to 

N c „ = | A„(a, + a?) + A n (a y + a*) ]©„ (73) 

N„ = [ A l2 (a, + a?) + A 22 (a, + <)]©„ (74) 

where the subscripted A -terms are elements of the 
inverse of the matrix in the left-hand side of matrix Eq. 
(72). Similarly, the expressions for the compliance 
coefficients Cj and C 2 that are needed for a buckling 
analysis and given by Eqs. (41)and (42), respectively, 
reduce to 


where the superscript (5) that identifies the central plate 
has been dropped for convenience and consistency with 
the previously used notation presented herein. The 
desired form of these two equations is obtained by 
expressing the temperature change in the restraining 
plates in terms of the temperature change in the central 

plate by 0r> k) = T (k) 0o (k = 1, 2, 3, and 4), and by defin- 
ing 


D oO) /|\ o(2) .a. 

afi = \ a Vi' + \-aft’ 

(67) 

k < 3 ) u< 4 > 

ah = V a $ + V 3 ® 

b b 

(68) 

i? = + ^a®x (21 

(69) 

= ~al 3) x ai + ^--a‘ 4) T (4) 
y b y b y 

(70) 


These expressions indicate that, in general, each plate 
can have a different uniform, isolated temperature field. 
When T (k) = 0 for all the restraining plates (k = 1, 2, 3, 
and 4), the restraining medium does not undergo any 
thermal deformations. 

The desired form of Eqs. (65) and (66) is given in 
matrix form by 



Next, because all the subcritical loads are zero-valued, 
N x (x,y) = -N5i, N y (x,y) = - N yl , and N, y (x,y) = N >yl = 0 
(see Fig. 2b). Substitution of these expressions into Eq. 
(71) yields 


(a,, + a") a„ 

^12 t ^22 ‘^■ 22 ) 



(72) 


which is a degenerate form of Eq. (31). Equation (72) 
shows that the thermally induced loads depend on the 


A n a, + A 12 q y 

A„(ot I + cC) + A I2 (a y + a*) 
A )2 a x + A 22 a y 

A t2 (a x + a, ) + A 22 (a y + a*) 


(75) 

(76) 


When the restraining medium is totally rigid, the com- 
pliances of the restraining medium in Eq. (72) are zero 

valued and the subscripted A -terms reduce to the usual 
membrane stiffnesses of laminated plate theory, which 
are denoted herein by subscripted A-terms. 

Examination of Eq. (72) indicates, as expected, that 
the Poisson-like lateral expansion or contraction of the 
restraining plates shown in Fig. 3 are totally indepen- 
dent; that is, this behavior is manifested by the absence 

ofa? 2 terms in Eq. (72). For a homogeneous, orthotropic 
restraining medium similar to that depicted in Fig. 2, this 
interaction of Poisson effects would be present, along 
with self-equilibrated shearing stresses and slight non- 
uniformities in the N x (x,y) and N y (x,y) distributions 
within the central plate. For this type of restraining me- 
dium, a reasonable approximation to the thermally in- 
duced loads is given by 


(a,, + a") (a 12 + a*) 1/ K,\ = / «, + a* \ 
(a 12 + a? 2 ) (a 22 + a 22 ) \ N y ,/ ^a y + a 


(77) 


where the compliances of the restraining medium could 
be estimated by subdividing the restraining medium into 
plates, neglecting the self-equilibrated shearing stresses, 
and applying a procedure similar to that used to deter- 
mine the compliances for the structural arrangement 
shown in Fig. 3. 

For the results that are presented subsequently in 
the present paper, it is convenient to introduce the rela- 
tive-compliance ratios R|=a?i/an and R 2 = a;n/a 22 * 
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Similarly, it is convenient to introduce the ratios of the 
coefficients of thermal expansion that are given by 

Ti = 0Cx/0tx and T 2 = oty/a y . Once a particular laminate 
is selected for investigation, the relevant properties of the 
restraining medium are specified in a relative manner, 
with respect to the central plate, by defining values for 

R,, R 2 , T,, and T 2 . With these ratios, the subscripted A- 
terms, appearing in Eqs. (73)-(76), become 

A„(l +R 2 )(A„A; 2 - A^) 

“ (1 +R,)(I + R J )A ll A a -Af 1 

t A| 2 (A M A n — A l2 ) 

,2 ~(1 + R,)( 1 +R 2 )A„A 22 - A^j 

A 22 (l +R,)(A„A 22 - a, 2 ) 

22 (1 +R,)(1 +R 2 )A„A 22 - A{ 2 

In addition, the compliance coefficients given by Eqs. 
(75) and (76) become 

1 + C, = 

(A„q, + A, ; q,)[(l +R,)(1 + R,)A„A„- A 2 ,,) 

|A„q,(l +R,)(1 + T,) + A 12 q,( 1 +T,))(A ll A n - A ! , : ) (8I) 


(78) 

(79) 

(80) 


1 +C 2 = 

(A„q, + A ;; q y )[ ( I + R,)(1 + Ri)A„A h - A,,] 

[A 12 q,( 1 +T,)+A :! q,(l +R,)(1 +T ; )|(A„A„- A 2 12 ) (82) 


Moreover, 

1 + C, = 

1 + C 2 

tA n a, + A l2 a y \|A l2 a,(l + T,) + A 22 a y (l + R,)(l + T 2 ) 1 
l A l2 a, + A 22 a y )[ A n a,( 1 + R,)(1 + T,)+A 12 a y (l +fjJ 


Results for Common Laminates and Discussion 


Results are presented in Figs. 4-17 that illustrate the 
effects of the inplane elastic resistance provided by the 
laboratory-type apparatus, previously described herein 
(see Fig. 3), on the restraint characteristics and buckling 
behavior of several common laminated plates subjected 
to a uniform temperature field. In particular, results are 
presented for several common balanced, symmetric lam- 
inates that are made of IM7/5260 graphite-bismaleimide 
material (see Table 1 ). The laminates include [(±45/0/ 
90)J S quasi-isotropic laminates, [(±45/0 2 )J s and [(±45/ 
90 2 ) J s quasi-orthotropic laminates, and [(±0)J S angle- 
ply laminates, where a positive value of the lamina fiber 
angle 0 is shown in Fig. 10. The [(±45/0 2 ) m ] s and [(±45/ 


90 2 ) m ] s laminates are described herein as quasi-orthotro- 
pic because of the presence of relatively small amounts 
of flexural anisotropy. All of the results are based on 
classical laminated-plate theory and the nominal ply 
thickness used in the calculations is 0.005 in. 

The buckling behavior of the [(±45/0/90)J s , 
[(±45/90 2 )J s , [(±45/0 2 ) J s , and [(±0) J s laminates is gov- 
erned by cases 1 through 4, that are defined by Eqs. (47) 
through (56), because balanced, symmetric laminates 
that are restrained against inplane movement do not ex- 
hibit shear loads when subjected to a uniform tempera- 
ture field. When a,/(l +Ci) is nonzero (see cases 1 
and 2), the key parameters for these laminates that are 

needed to determine their stability are &|/(1 + Ci) and 

( 1 + Cl , or equivalently, 60/(1 +Ci) and 
\l + C2JG1 

& 2 /(l +C 2 ). When dj/(l + Ci) = 0, the key parameter is 
ft 2 /( 1 + C 2 ), as defined by cases 3 and 4. 

Results are presented in Figs. 4-6 that show the 
nondimensional compliance coefficient (1 +C 1 )(1 + T,) 
for the [(±45/0/90)J s , [(±45/90 2 ) J s , and [(±45/0 2 )J s bal- 
anced, symmetric laminates, respectively. Specifically, 
results are shown as a function of the relative-compli- 
ance ratio R t = a* /an, for discrete values of the relative- 

compliance ratio R 2 = a 22 /a 22 that range from 0 to 0.5. 
For these results, the ratios of the coefficients of thermal 

expansion, that are given by Ti=(Xx/tx x and 

T 2 = ay/(Xy, are equal; that is T, = T 2 . Values of R, = R 2 
= 0 correspond to plate that is rigidly restrained against 
thermal expansion or contraction. Nonzero values of T, 
account for thermal expansion of the restraining medi- 
um. The parameter (1 + C, )( 1 + T,) is used in the figures 
because C, was found to be independent of the number 
of laminates plies (determined by the symbol m in the 
stacking sequence specification) and because Eq. (81) 
can be expressed in a manner that is independent of T t 
for the special case of Tj = T 2 ; that is, 

(1 +C,)(1 +T,) = 

(A„a, + A 12 a y )[(l +R ,)(1 +R 2 )A,,A 22 - A' ; ] 

[A n a,(l +Rd + A l 2 a y ](A n A 22 -A 2 12 ) (84) 

Thus, graphs of 6 l i, that are given in Ref. 7 for specific 
values of the stacking sequence number m, can be used 
with the graphs of (1 + C,)(l + T,) provided herein to 
determine the buckling behavior for a large number of 
restraint configurations that include thermal deforma- 
tions of the restaining plates. 


12 

American Institute of Aeronautics and Astronautics 


The results in Figs. 4-6 indicate a strong effect of 
laminate stacking sequence on the value of (1 + CjXl + 
Tj), with the axially stiff [(±45/0 2 ) m ] s laminates exhibit- 
ing the widest variation with variations in Rj and R 2 . Re- 
sults for the transversely stiff [(±45/90 2 ) m ] s laminates 
exhibit the least variation. Together, the results in Figs. 
4-6 indicate that the parameter (1 + C,)(l + Tj) is a pos- 
itive-valued, linear function of R! and increases with in- 
creases in R, or R 2 . Because (1 + Cj)(l + T { ) > 0 and 

6c i > 0 for these laminates (see Fig. 3 of Ref. 7), it fol- 
lows that 6tj/( 1 + C i) > 0 as long as (1 + T,) > 0. Thus, 
the buckling behavior of these laminates is governed by 
case 1 for uniform heating and by case 2 for uniform 
cooling when (1 + T,) > 0, which typically includes most 
situations of practical importance. 

Results are presented in Fig. 7 that show the nondi- 
mensional compliance coefficient ratio (1 +0^/(1 + C 2 ) 
as a function of the relative-compliance ratio R,, for dis- 
crete values of the relative-compliance ratio R 2 that 
range from 0 to 0.5. These results are also for the case 
where the ratios of the coefficients of thermal expansion 
are equal (Tj = T 2 ). Moreover, the results are also inde- 
pendent of the number of laminate plies and the value of 
T,. The black solid lines, the gray solid lines, and the 
gray dashed lines in the figure are for the [(±45/0/90) m ] s , 
[(±45/0 2 ) m ] s , and [(±45/90 2 )J s laminates, respectively. 
The expression for (1 + C^/0 + C 2 ) that is obtained by 
simplying Eq. (83) accordingly is given by 

i +c, f A n a, + A <; a y \ 

1 +C l “VA„a. + A*aJ 

Corresponding values of ft 2/61 that are needed to com- 
pute the parameter ( I + 1" are given in Fig. 4 of 

Ref. 7, for specific numbers of laminate plies, and are all 
positive valued. 

The results in Fig. 7 also show a strong effect of 
laminate stacking sequence on the value of (1 + C,)/(l + 
C 2 ). The greatest sensitivity to variations in R p as indi- 
cated by the slopes of the lines, is exhibited by the axially 
stiff [(±45/0 2 ) m ] s laminates, followed by the 
[(±45/0/90) m ] s quasi-isotropic laminates and then the 
transversely stiff [(±45/90 2 )J s laminates. This trend ap- 
pears to be reversed with respect to variations in R 2 . In 
addition, Eq. (85) and the results in Fig. 7 indicate that 
the ratio (1 + Cj)/(1 + C 2 ) is a positive-valued linear func- 
tion of R, for all three laminates. The values of (1 + Cj)/ 
(1 + C 2 ) for the axially stiff [(±45/0 2 )J s laminates in- 
crease with increases in R, or R 2 . In contrast, the values 


A n a t + A 22 a y ( 1 + R,) 


A,,a x ( 1 + R 2 ) + A l2 a, 


(85) 


of (I + C,)/(l + C 2 ) for the [(±45/0/90) J s quasi-isotropic 
laminates increase with increases in R, and decrease with 
increases in R 2 . Likewise, the values of (1 + C^/0 + C 2 ) 
for the transversely stiff [(±45/90 2 )J s laminates decrease 
with increases in Rj or R 2 . 

Results that are similar to those presented in Figs. 
4-7 are presented in Figs. 8 and 9 for the special case 
where the relative-compliance ratios are equal and the ra- 
tios of the coefficients of thermal expansion are equal; 
that is R, = R 2 and T, = T 2 . The black solid lines, the gray 
solid lines, and the gray dashed lines in the figures are for 
the [(±45/0/90)J s , [(±45/0 2 U, and [(±45/90 2 )J s lami- 
nates, respectively. Additionally, the results are inde- 
pendent of the number of laminate plies and the results in 
Fig. 9 are also independent of the value of T,. 

Like the results in Figs. 4-6, the results in Fig. 8 
also show a strong effect of laminate stacking sequence 
on the value of (1 + C,)(l + T,) for the axially stiff 
[(±45/0 2 ) m ] s laminates. The results for these axially stiff 
laminates exhibit the greatest variations in (1 + C,)(l + 
T,) with variations in R { and a slight nonlinearity in the 
curve. The results for the [(±45/0/90)J s quasi-isotropic 
laminates and the transversely stiff [(±45/90 2 ) m ] s lami- 
nates exhibit the least variation with R } and are nearly 
identical. For all three laminate types, the parameter (1 
+C,)(1 +Tj) is a positive-valued, monotonically increas- 
ing function of the relative-compliance ratio R,. 

The results in Fig. 9 also show a strong effect of 
laminate stacking sequence on the value of ( 1 + Cj)/(1 + 
C 2 ), like the results in Fig. 7. The greatest sensitivity to 
variations in R l is exhibited by the axially stiff 
[(±45/0 2 ) m ] s laminates, which exhibit a slightly nonlin- 
ear, monotonically increasing variation with R,. The 
next greatest variation in (1 + C j )/( 1 + C 2 ) with R t is ex- 
hibited by the transversely stiff [(±45/90 2 ) m ] s laminates, 
which exhibit a slightly nonlinear, monotonically de- 
creasing variation with R,. The least sensitivity to vari- 
ations in R, is exhibited by the [(±45/0/90) m ] s quasi- 
isotropic laminates, for which (1 + C,)/(1 +C 2 )= 1 for all 
values of R,. This constant value of unity is seen by ex- 
amining Eq. (85) and noting that A n = A 22 and a x = a y 
for this family of laminates. 

In Ref. 7, it is shown that the values of fti and 

ft 2 /ft 1 are all positive for the [(±45/0/90) m ] s , [(±45/0 2 ) J s , 
and [(±45/90 2 ) m ] s laminates. Likewise, the results in 
Figs. 4-9 indicate that (1 + Cj)/(1 + C 2 ) > 0 and that (1 + 
C,) > 0 as long as (1 + T,) > 0. Thus, for this practical 
range of T p these laminates are loaded by biaxial com- 
pression when uniformly heated and by biaxial tension 
when uniformly cooled. Therefore, these laminates can 
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buckle only when subjected to uniform heating (see case 
1) when (1 + T,) > 0. 

Results are presented in Figs. 10 and 1 1 that show 
the stiffness-weighted thermal expansion parameters, 

d+cnn+ T .) and (TTclW respectively, for the 
[(±0)J S angle-ply laminates. The results in these figures 
are presented as a function of the fiber angle 0 and are 
for the special case where the relative-compliance ratios 
are equal and the ratios of the coefficients of thermal ex- 
pansion are equal (R, = R 2 and T, = T 2 ). The black solid 
lines correspond to different values of the relative-com- 
pliance ratio R,, that range from values of 0 to 0.5. Ad- 
ditionally, the results are independent of the number of 

6 

laminate plies. An expression for ^ +ci)(i +Ti ~ ) ls °^“ 
tained by combining Eqs. (38) and (81) and then simply- 
ing. The result is given by 

6t, __ t^AnCX.O +R t ) + A l2 q y ](A„A 22 - A 2 2 ) /o ^ 

(l+C.JJI+T,) 12yD^[(l + R,) 2 A u A 22 -A 2 l2 ] 

Similarly, combining Eqs. (39) and (82) and then sim- 
plying gives 

t 2 [A 12 q t + A 22 q y (i R ,)]( A n A 23 - A, 2 ) 

(1 + C 2 |(1 +T,) 12D 22 [(l + Ri) 2 A,,A I2 - A 2 l2 ] 

The curves in Figs. 10 and 1 1 show a very strong 
dependence on the fiber angle 0. In particular, 

n — rr-xTi — =— r is shown to increase, for the most part, with 
(1 +Ci)(l +Ti) 

ft 

increasing values of 0, and ^ +c 2 ~ )( 2 i T ti) ]S s ^own to de- 
crease. In addition, both stiffness-weighted thermal ex- 
pansion parameters are shown to decrease with 
increasing values of the relative-compliance ratio Rp 
For the curves that correspond to Rj = 0, 0.1, 0.2, and 
0.3; both stiffness- weighted thermal expansion parame- 
ters are positive for all values of the fiber angle 0, which 

means that &i/(l + Ci)> 0 as long as (1 + T { ) > 0. For 
these values of R p and for (1 + T,) > 0, the laminates are 
loaded by biaxial compression for uniform heating (case 
1) and biaxial tension for uniform cooling (case 2). 
Thus, buckling can occur only for uniform heating, as 
governed by case 1 . 

The results for the angle-ply laminates with R, = 
0.4 and 0.5 that are shown in Figs. 10 and 1 1 indicate the 
possibility of several thermally induced loading condi- 
tions, depending on the particular value of the fiber angle 
0. First, the curves that correspond to R l = 0.4 and 0.5, 

have zero values for jz — that are indicated in 

{1 +C|)(1 +T t ) 


Fig. 10, by the unfilled circular symbols. Specifically, 
the curve that corresponds to R ] = 0.4 has zero values at 
approximately 0 = 18 deg and 29 deg. The curve that 
corresponds to R x = 0.5 has zero values at approximately 
0=15 deg and 32 deg. For these specific laminate and 

restraint configurations with p +c^){\ + Fig. 11 in- 

dicates that the corresponding values of ^ + c ^ are 
positive as long as (1 + Tj) > 0, which means that the 
plates are loaded by only transverse compression or ten- 
sion for uniform heating or uniform cooling, respective- 
ly. Thus, the plates may buckle only under uniform 
heating, into a wide-column mode, as governed by Eq, 
(54b) of case 3. Second, the curve that corresponds to Ri 

= 0.4 has positive values for ^ +c ^ 1 1 and 

(l +€ 2)0 +f V ) ^ or va ^ ues (approximately) 0 < 18 deg, 
29 < 0 < 61 deg, and 0 > 72 deg as long as ( 1 + Tj) > 0. 

Third - ( , + c i)('i + t i ~ ) is ne 8 ative and (i + c ^( 2 | + T|) * S P° si - 

live for 1 8 < 0 < 29 deg, and (1+c ^ +T]) is positive and 

(i + 02)0 +T|j * s ne S at * ve f° r 61 < 0 < 72 deg as long as (1 
+ Tj) > 0. Thus, for the range 0 < 18 deg, 29 < 0 < 61 
deg, and 0 > 72 deg; the laminates are loaded by biaxial 
compression for uniform heating (case 1 ) and biaxial ten- 
sion for uniform cooling (case 2) and, thus, buckling can 
occur only for uniform heating. For 1 8 < 0 < 29 deg, the 
laminates are loaded by axial tension and transverse 
compression for uniform heating (case 2) and oppositely 
directed loads for uniform cooling (case 1). Thus, for 1 8 
< 0 < 29 deg, buckling can occur for uniform heating or 
cooling. For 61 < 0 < 72 deg, the laminates are loaded 
by axial compression and transverse tension for uniform 
heating (case 1) and oppositely directed loads for uni- 
form cooling (case 2). Thus, for 61 < 0 < 72 deg, buck- 
ling can also occur for uniform heating or cooling. 

The results in Figs. 10 and 1 1 for the [(±0) m ] s an- 
gle-ply laminates with R, = 0.5 and (1 + T,) > 0 indicate 
behavioral trends that are similar to those previously de- 
scribed for the laminates with R t = 0.4. In particular, for 
the range 0 < 15 deg, 32 < 0 < 58 deg, and 0 > 75 deg; 
the laminates are loaded by biaxial compression for uni- 
form heating (case 1) and biaxial tension for uniform 
cooling (case 2), and buckling can occur only for uni- 
form heating. For 15 < 0 < 32 deg, the laminates are 
loaded by axial tension and transverse compression for 
uniform heating (case 2) and oppositely directed loads 
for uniform cooling (case 1). Thus, for 15 < 0 < 32 deg, 
buckling can occur for uniform heating or cooling. For 
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58 < 0 < 75 deg, the laminates are loaded by axial com- 
pression and transverse tension for uniform heating (case 
1) and oppositely directed loads for uniform cooling 
(case 2). Thus, for 58 < 0 < 75 deg, buckling can occur 
for uniform heating or cooling. 

Values for the critical temperature changes for the 
[(±45/0/90)J s , [(±45/0 2 )J s , and [(±45/90 2 )J s laminates 
and for the [(±0)J S laminates with m = 1 and m > 5 are 
shown in Figs. 12-17. In these Figures, the buckling re- 
sults are based on the solutions for infinitely long plates 
and uniform heating. The solid and dashed curves corre- 
spond to results for plates with simply supported and 
clamped flexural boundary conditions, respectively. In 

particular, Fig. 1 2 shows values of for the [(±45/ 

0/90)J s , [(±45/0 2 )J s , and [(±45/90 2 )J s laminates, that 
are fully (rigidly) restrained against thermal expansion or 
contraction, as a function of the number of laminate 
plies. This fully restrained condition is indicated by the 

symbol 0o\ Similarly, Fig. 13 shows values of 
for the [(±0) m ] s laminates, with m = 1 and m > 5, that 
are fully (rigidly) restrained against thermal expansion or 
contraction. The results in Figs. 14-16 show the critical- 


temperature-change ratio as a function of the 

0o(t + C i ) 

nondimensional compliance ratio (1 + C t )/( 1 + C 2 ) and 
the stacking sequence numbers m = 1 and m = 8 for the 
[(±45/0/90)J s , [(±45/0 2 )J s , and [(±45/90 2 )J s laminates, 
respectively. Likewise, the results in Fig. 17 show the 


e8 r 


as a func- 


critical-temperature-change ratio 

00 (1 + Ci) 

tion of the nondimensional compliance ratio (1 + Cj)/(1 
+ C 2 )forthe [(±0) m ] s laminates, with m > 5, and for val- 
ues of the fiber angle 0 = 30, 45, and 60 deg. In Figs. 
14-17, the results are based on the presumption that (1 + 


C() > 0, such that ft i0o/(l + Ci) > 0 and the buckling 
responses are governed by case 1 (Lj = 1 and L 3 = 0) for 
uniform heating. This presumption is reasonable be- 


cause fti and ft/ft, are all positive for these laminates, 

which means that [ * + > 0 for the range of (1 + C,)/ 

V l + C2/G1 

(1 + C 2 ) selected for the independent variable. More- 
over, the parameter (1 + C,) appears in the expression 


Qft r 

©o(l +Ci) 


as a value to be selected. 


The results in Figs. 14-17 give the critical tempera- 
ture change for a wide range of restraint conditions, rel- 
ative to the critical temperature change for the 
corresponding plate that is fully (rigidly) restrained 
against thermal expansion or contraction (Figs. 12 and 
13). For example, values of (t + C,) and (1 + C,)/(l + 


C 2 ) for the [(±45/0/90) m ] s quasi-isotropic laminates can 
be obtained from Figs. 4 and 7, respectively, and used 
with Fig. 14 to obtain critical temperature changes for a 
very large number of restraint configurations that are de- 
fined by the relative-compliance parameters Rj and R 2 
and by the ratio of the coefficients of thermal expansion 
given by T r This approach was used in the present study 
because the results for a given fully restrained plate serve 
as a reasonable reference point for assessing the effects 
of a given elastic restraining medium. 

The buckling results shown in Figs. 12 and 13 for 
the fully restrained laminates were obtained by comput- 
ing the ratios ft 2 /ft 1 first, in accordance with the proce- 
dure of case 1. Next, the load factor L 2 was determined 
by using Eq. (47) with C ] = 0 and C 2 = 0. The load 
factor L 3 is zero-valued, because balanced laminates do 
not induce shear loads when subjected to restrained ther- 
mal expansion or contraction. Then, the critical value of 

the loading parameter p cr was determined for each of the 
laminates and the critical temperature change was ob- 
tained by using Eq. (49) with C, = 0. Equations (47) and 

(49) were used because fti0 () >O for these laminates 
when subjected to uniform heating. For these laminates, 
the critical value of the loading parameter depends on the 
flexural boundary conditions, the plate flexural orthotro- 
py and flexural anisotropy, and the value of the load fac- 
tor L 2 , which depends on the stiffness-weighted laminate 
thermal-expansion parameters. The stiffness-weighted 
laminate thermal-expansion parameters depend on the 
plate membrane orthotropy and the plate stiffnesses as- 
sociated with pure bending action. 

For the buckling results in Figs. 14-17, the nondi- 
mensional compliance ratio (1 + Cj)/( I + C 2 ) was treated 
as an independent variable. For a fixed value of this ra- 
tio, buckling results were obtained by again computing 

the ratios ft/ft, first and then the load factor L 2 by using 
Eq. (47), in accordance with the procedure of case 1. 
Then, the critical value of the loading parameter p cr was 
determined from a buckling analysis for each of the lam- 
inates and the critical temperature change was obtained 
by using Eq. (49) in the form 


12b ; Q 0 _ Per 

7tV(l+C,) 


( 88 ) 


where per = pcr(P, Y, 8, L 2 ) for a given sel of flexural 
boundary conditions and the load factor L 2 depends on 
the stiffness-weighted laminate thermal-expansion 
parameters and the nondimensional compliance ratio 
(1 + Cj)/(1 + C 2 ). Next, results obtained from Eq. (88) 
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were divided by the corresponding results for the corre- 
sponding laminates that are fully restrained (Cj ~ C 2 - 
0), which are given symbolically by 


12bj5>«_P« 
! — 


nv 


a, 


(89) 


Here, it is reiterated that the fully restrained condition is 
indicated by the symbol 0". Similarly, the symbol p" r 
is used to denote the critical value of the loading param- 
eter for a given fully restrained laminate. For this situa- 
tion, per = pcr(P, Y, 8, L”), where L" = ft 2 /Gi • Finally, 
dividing Eqs. (88) and (89) gives the result 


e;; _ % 

©:(i+c,) p: 


(90) 


It is useful to point out, that when p cr and p? r corre- 
spond to wide-column buckling modes, Eq. (90) reduces 
to 

&: n+c,)-' 

©“( 1 + C,) v 1 +C 2 ) (9D 


This special case is indicated in Figs. 14-17 by the solid 
gray curve in each figure (the solid gray curve lies 
underneath the solid black curve in Fig. 14). 

The results shown in Fig. 12 for the fully restrained 
plates indicate that the critical temperature change 

^y0o is highly dependent on the arrangement of the 

45-deg, 0-deg, and 90-deg plies. As the number of plies 
increases, the magnitude of the critical temperature 
change decreases for the clamped and simply supported 
[(±45/0 2 ) m ] s laminates and increases for the clamped and 
simply supported [(±45/0/90)J s laminates and the sim- 
ply supported [(±45/90 2 )J s laminates. The critical tem- 
perature change varies only slightly for the clamped 
[(±45/90 2 ) m ] s laminates. Overall, the plates with 
clamped edges are more buckling resistant than the cor- 
responding plates with simply supported edges, as ex- 
pected. The clamped and simply supported [(±45/90 2 )J s 
laminates require the most heating to cause buckling and 
the simply supported [(±45/0 2 )J s laminates require the 
least amount of heating. 

The results shown in Fig. 13 for the fully restrained 
plates indicate that the critical temperature change 

^©o for the [±0] s and [(±0)J S (m > 5) laminates is 

highly dependent on the fiber angle 0. In particular, the 
results for the clamped and simply supported [±0] s and 
[(±0) m ] s (m > 5) laminates show substantial increase in 


the critical temperature change with increasing 0 for val- 
ues up to approximately 55 deg and 62 deg, respectively, 
followed by a significant decrease. In addition, the 
[(±0)J S (m > 5) laminates require more heating to cause 
buckling than the corresponding [±0] s laminates. Over- 
all, the plates with clamped edges are much more buck- 
ling resistant than the corresponding plates with simply 
supported edges, as expected. 

The data presented in Figs. 1 2 and 1 3 are given in a 
compact form and are applicable to an infinite range of 
plate width-to-thickness ratios, b/t. Once the number of 
laminate plies is selected, the critical temperature change 

0o can be found as a function of the plate width b for 
a given laminate family. However, it is important to 
keep in mind the limitations of classical laminated -pi ate 
bending theory as the plate width-to-thickness ratio b/t 
becomes smaller than a value of approximately 20. Re- 
sults of this type are useful in structural design and are il- 
lustrated for the fully restrained [(±45/0/90)J s quasi- 
isotropic laminates in Fig. 13 of Ref. 7. 

The results in Figs. 14-17 indicate a general, com- 
mon trend of a reduction in the critical-temperature- 

change ratio — - with increases in the nondimen- 

sional compliance ratio (1 + C j )/ ( 1 + C 2 ), for all the lam- 
inates considered and both sets of flexural boundary 
conditions. This reduction in (1 + C j )/( 1 + C 2 ) corre- 
sponds to a reduction in the thermally induced transverse 
compression load and hence the load factor L 2 . For many 
of the laminate configurations shown in the figures, the 
critical-temperature-change ratios are coincident with 
points of the gray curves that are defined by Eq. (9 1 ). In 
Fig. 14, the solid black curve shown is on top of the gray 
curve defined by Eq. (91). This coincidence indicates 
that an infinitely long laminate that is either fully re- 
strained or elastically restrained exhibits a wide-column 
buckling mode for both restraint cases. In contrast, the 
curves shown in Figs. 14-17 that are not coincident with 
the solid gray curve, which are mostly for clamped 
plates, indicate the absence of a wide-column buckling 
mode. In addition, the effect of varying the number of 
plies in the [(±45/0/90)J s , [(±45/0 2 )J s , and [(±45/90 2 ) m ] s 
laminates varies significantly with the laminate type. 
For plates with finite length, these trends are expected to 
change dramatically because moderate-aspect-ratio, fi- 
nite-length plates do not exhibit wide-column buckling 
modes. 


Generic Results and Examples 

In Ref. 7, generic results are presented that enable 
one to obtain critical temperatures for a wide range of 
laminate configurations that are subjected to fully re- 
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strained thermal expansion or contraction by using 
known buckling-behavior results for anisotropic plates 
subjected to combined mechanical loads. The work pre- 
sented in Ref. 7 is based on the simplicity of equations 
like Eq. (49) and consists of data that define fundamental 
parameters for a wide variety of selected laminate con- 
structions. Specifically, figures that show the stiffness- 

weighted laminate thermal-expansion parameters <Xi, 

(x 2 , and <5t3 ; the flexural orthotropy parameter (3; and the 
flexural anisotropy parameters y and 8 as a function of 
material system and laminate stacking sequence are pre- 
sented in Ref 7. In addition, figures that show the criti- 
cal value of the loading parameter p cr , or buckling 
interaction curves, as a function of flexural boundary 
conditions, (3, y, 8, L,, L 2 , and L 3 are also presented for a 
broad range of parameters and combined mechanical 
loading conditions. Results of this type are significant 
because they illustrate the key aspects of the behavior 
and show overall trends and sensitivity of the behavior to 
changes in the parameters. Moreover, these results en- 
able one to apply knowledge of, and experience with, the 
behavior of anisotropic plates subjected to mechanical 
loads to the problem of restrained thermal expansion or 
contraction. 

The present study builds upon the work in Ref. 7 by 
presenting generic data in this section that can be used 
with the results in Ref. 7 to determine the effect of an 
elastic restraining medium, like the one shown in Fig. 3, 
on the thermal buckling behavior for a very wide range 
of laminate constructions and restraint configurations. 
Specifically, results are presented that show the effects 
of material systems and laminate stacking sequence on 
the fundamental parameters that are used to determine 
the thermally induced mechanical loads. Moreover, cor- 
responding values of the fundamental nondimensional 
parameters (3, y, and 8 that are required to use the generic 
mechanical-buckling-behavior results that appear in Ref. 
7 are also presented. In addition, the procedure used to 
determine the critical temperature change by using the 
generic results is discussed. 

Values of the Fundamental Parameters 

Values of the fundamental parameters that are used 
herein are presented in this section for transversely stiff 
[(±45/90 2 )J s laminates, [(±45/0/90)J s quasi-isotropic 
laminates, axially stiff [(±45/0 2 )J s laminates, and 
[(±0)J S angle-ply laminates. Nine different contempo- 
rary material systems are used. These material systems 
include boron-aluminum, S-glass-epoxy, a typical bo- 
ron-epoxy, AS4/3501-6 graphite-epoxy, AS4/3502 
graphite-epoxy, IM7/5260 graphite-bismaleimide, Kev- 
lar 49-epoxy, IM7/PETI-5, and P-100/3502 pitch-epoxy 
materials. The mechanical properties of these material 


systems are presented in Table 1 and the nominal ply 
thickness is 0.005 in. 

Parameters for [(±45/9Q 2 )J s laminates Values of 
(1 +Cj)(l + T,) and (1 +C,)/(1 + C 2 ) for [(±45/90 2 )J s 
laminates are presented in Figs. 18 and 19, respectively, 
as a function of the relative-compliance ratio Rj for the 
nine material systems. The graph in Fig. 19b is an en- 
largement of the portion of the graph in Fig. 19a that con- 
tains the closely grouped curves. The results in these 
figures are for the special case where the relative-compli- 
ance ratios are equal and the ratios of the coefficients ot 
thermal expansion are equal (Rj = R 2 and T ] = T 2 ), and 
are independent of the number of laminate plies. The re- 
sults show a wide variation in (1 + C,)(l + T t ) and (I + 
C , )/( 1 + C 2 ) with material system. Moreover, (1 + C,)(l 
+ T,) is positive for all the material systems, and the 
curves shown in Fig. 18 increase monotonically with in- 
creases in the relative-compliance ratio R,, except the 
curve for the P-100/3502 pitch-epoxy material. Similar- 
ly, (1 + C,)/(l + C 2 ) is positive for all the material sys- 
tems and the curves shown in Figs. 1 9a and 19b decrease 
monotonically with increases in the relative-compliance 
ratio R, , except the curve for the Kevlar 49-epoxy mate- 
rial. In addition, the magnitudes of (1 + C , )/( 1 + C 2 ) for 
the Kevlar 49-epoxy material are much greater than 
those for the other eight materials. The curves for the bo- 
ron-aluminum, S-glass-epoxy, and boron-epoxy materi- 
als exhibit values of (1 + C,)/(l + C 2 ) that are nearly 
equal to unity, which indicates a relatively small contri- 
bution of the elastic restraining medium to the thermally 
induced transverse load (y-coordinate direction shown in 
Fig. 2b). 

Values of and & 2 /&\ for [(±45/90 2 )J s lami- 
nates are presented in Figs. 20 and 21, respectively, as a 
function of the number of laminate plies for the nine ma- 
terial systems. The graph in Fig. 21b is an enlargement 
of the portion of the graph in Fig. 21a that contains the 
closely grouped curves. The results show a wide varia- 
tion in with material system, a relatively benign vari- 
ation with the number of plies, and that (x \ is positive for 

all the material systems. However, is nearly zero val- 
ued for the P-100/3502 pitch-epoxy material. Likewise, 
the results in Figs. 21a and 21b show a substantial varia- 
tion in 6i 2 /0i\ with material system, and that & 2 /&i is 
positive for all the material systems except the P-100/ 
3502 pitch-epoxy and Kevlar 49-epoxy materials. Spe- 
cifically, the values for the Kevlar 49-epoxy material are 
slightly negative and nearly equal to zero and those for 
the P-100/3502 pitch-epoxy material are negative and 
have much larger magnitudes than the corresponding 
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values for the other materials. 

The results in Figs. 18 and 20 indicate that the val- 
ues of (1+c ^ 1 1+T - -y are positive for the [(±45/90 2 )J s 
laminates, for all nine material systems, which means 
that i/(l + C i) > 0 as long as (1 + T,) > 0. Similarly, 
the results in Figs. 19 and 21 indicate that the values of 

( 1 + C 1 arC ^ or mater ^ a ^ s exce pt the P- 

100/3502 pitch-epoxy and Kevlar 49-epoxy materials. 
Thus, for the majority of the materials, the laminates are 
loaded by biaxial compression when uniformly heated 
and by biaxial tension when uniformly cooled, and as a 
result can buckle only for uniform heating as long as (1 
+ T,) > 0 (see case 1 ). For the Kevlar 49-epoxy and P- 
100/3502 pitch-epoxy materials, the laminates are load- 
ed by axial compression and transverse tension when 
uniformly heated and by oppositely directed loads when 
uniformly cooled, as long as (1 + Tj) > 0. For this situa- 
tion, the laminates can buckle for uniform heating (case 
1) or uniform cooling (case 2). Overall, Figs. 18-21 
show that there is a wide range of possibilities available 
for tailoring the thermal buckling characteristics of a 
laminate family and restraining medium. 

The remaining parameters needed for a buckling 
analysis of the long transversely stiff [(±45/90 2 )J s lami- 
nates are the the nondimensional flexural orthotropy pa- 
rameter P and the flexural anisotropy parameters y and 
8. Values for these parameters are presented in Figs. 22- 
24 as a function of the number of laminate plies, for the 
nine material systems. The results in Fig. 22 show a se- 
ries of monotonically decreasing curves that each ap- 
proach a constant value of P from above as the number 
of plies increases. Similarly, the results in Figs. 23 and 
24 show a series of monotonically decreasing curves that 
approach negligible values for y and 5 , respectively, 
from above as the number of plies increases. The largest 
and smallest values of p, y and 8 are exhibited by the 
laminates made of the P-100/3502 pitch-epoxy and bo- 
ron-aluminum materials, respectively. 

Parameters for [(±45/0/90)J s laminates Values of 
(1 + Cj)(l + T,) for [(±45/0/90)J s quasi-isotropic lami- 
nates are presented in Fig. 25 as a function of the rela- 
tive-compliance ratio Rj for the nine material systems. 
The results in this figure are also for the special case 
where the relative-compliance ratios are equal and the ra- 
tios of the coefficients of thermal expansion are equal (Rj 
= R 2 and Tj = T 2 ), and are independent of the number of 
laminate plies. Unlike the results for the [(±45/90 2 ) m ] s 
laminates, the results in Fig. 25 show a relatively small 
variation in ( 1 + C,)( 1 + T,) with material system for the 
quasi-isotropic laminates. Moreover, (1 + C,)(l + Tj) is 


positive for all the material systems, and the curves 
shown in Fig. 25 increase monotonically with increases 
in the relative-compliance ratio R,, with no exceptions. 
Values for (1 + Cj)/(1 + C 2 ) were determined to be equal 
to unity, for all nine material systems and for the values 
of the relative-compliance ratio R, shown in Fig. 25. 
This result indicates that there is no contribution of the 
elastic restraining medium to the thermally induced 
transverse load (y-coordinate direction shown in Fig. 
2b). 

Values of &\ and (iiKx \ for [(±45/0/90) m ] s quasi- 
isotropic laminates are presented for the nine material 
systems in Figs. 14 and 15, respectively, in Ref. 7. The 

results also show a wide variation in frj with material 

system for the quasi-isotropic laminates. Moreover, 
is positive for all the material systems except the P-100/ 
3502 pitch-epoxy material. Similarly, & 2/&1 is positive 
for all the material systems. 

The results in Fig. 14, in Ref. 7, and Fig. 25 herein 
G 

indicate that the values of +Tl ) are P ositive for 

the [(±45/0/90)J s quasi-isotropic laminates, for all nine 
materials, except the P-100/3502 pitch-epoxy material. 

Thus, <5c i/( 1 + C 1 ) > 0 for alt materials but the P-100/ 
3502 pitch-epoxy material as long as (1 + T,) > 0. Sim- 
ilarly, the results in Fig. 15, in Ref. 7, and the result (1 + 

C , )/( 1 + C 2 ) = 1 indicate that the values of | j ^ 2 )^ are 

positive for all the materials. Thus, for all materials ex- 
cept the P-100/3502 pitch-epoxy material, the [(±45/0/ 
90) m ] s quasi-isotropic laminates are loaded by biaxial 
compression when uniformly heated and by biaxial ten- 
sion when uniformly cooled, and as a result can buckle 
only for uniform heating as long as (1 + Tj) > 0 (see case 
1). For the P-100/3502 pitch-epoxy material, the lami- 
nates are loaded by biaxial tension when uniformly heat- 
ed and by biaxial compression when uniformly cooled. 
These laminates must be cooled to buckle (case 1 ). 

Values of the nondimensional orthotropy parame- 
ter p and the nondimensional anisotropy parameters y 
and 8 needed for a buckling analysis of the [(±45/0/ 
90) m ] s quasi-isotropic laminates made of the same 9 ma- 
terial systems have been presented graphically in Ref. 6. 
The parameters are presented graphically in Ref. 6 as a 
function of the stacking-sequence number m, along with 
a discussion of their characteristics. The results in Ref. 
6 for these parameters show a series of monotonically 
decreasing curves that approach P = 1, y = 0, and 8-0 
from above as the number of plies increases. A homoge- 
neous, isotropic material has values of p = 1, y= 0, and 
8 = 0. Thus, the results for P, y, and 8 in Ref. 6 give, to 
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some extent, a quantitative measure of quasi-isotropy. 

Parameters for [(±45/Q 2 )J s laminates Values of 
(1 +Cj)(l + Tj) and (1 + C,)/(1 + C 2 ) for axially stiff 
[(±45/0 2 )J s laminates are presented in Figs. 26 and 27, 
respectively, as a function of the relative-compliance ra- 
tio R, for the nine material systems. The results in these 
figures are also for the special case where the relative- 
compliance ratios are equal and the ratios of the coeffi- 
cients of thermal expansion are equal (R x = R 2 and T x = 
T 2 ). Moreover, the results are independent of the number 
of laminate plies. Like for the [(±45/90 2 ) m ] s laminates, 
the results in Figs. 26 and 27 for the [(±45/0 2 ) m ] s lami- 
nates show a wide variation in (1 + C,)(l + Tj) and (1 + 
C,)/(l + C 2 ) with material system. Moreover, (1 + 0^(1 
+ T,) is positive for all the material systems, and the 
curves shown in Fig. 26 increase monotonically with in- 
creases in the relative-compliance ratio R p except the 
curve for the Kevlar 49-epoxy material. Similarly, (1 + 
Cj)/(1 + C 2 ) is positive for all the material systems and 
the curves shown in Fig. 27 increase monotonically with 
increases in the relative-compliance ratio R { , except the 
curve for the Kevlar 49-epoxy material. In addition, the 
magnitudes of ( 1 + Cj)/( 1 + C 2 ) for the Kevlar 49-epoxy 
material are generally much smaller than those for the 
other eight materials. The trends exhibited by the curves 
in Fig. 27 for the [(±45/0 2 )J s laminates are generally op- 
posite of the trends exhibited in Figs. 1 9a and 1 9b for the 
corresponding [(±45/90 2 )J s laminates. The curves 
shown in Fig. 27 for the boron-aluminum, S-glass-ep- 
oxy, and boron-epoxy materials exhibit values of (1 + 
C,)/(l + C 2 ) that are nearly equal to unity, which indi- 
cates a relatively small contribution of the elastic re- 
straining medium to the thermally induced transverse 
load (y-coordinate direction shown in Fig. 2b). 

Values of ftt and <5 c 2 /Gi for [(±45/0 2 )J s laminates 
are presented in Figs. 28 and 29, respectively, as a func- 
tion of the number of laminate plies for the nine material 
systems. The graph in Fig. 29b is an enlargement of the 
portion of the graph in Fig. 29a that contains the closely 
grouped curves. The results for these axially stiff lami- 
nates also show a wide variation in ft i with material sys- 
tem, and a generally benign variation with the number of 

plies, and that fti is positive for all the materials except 
the P-100/3502 pitch-epoxy and Kevlar 49-epoxy mate- 
rials. Moreover, fti is nearly zero valued for the Kevlar 
49-epoxy material. Likewise, the results in Figs. 29a and 
29b show a substantial variation in ft 2 /fti with material 

system, and that ft 2 /fti is positive for all the material 
systems except the P-100/3502 pitch-epoxy and Kevlar 
49-epoxy materials. Specifically, the values for the P- 


100/3502 pitch-epoxy material are slightly negative and 
nearly equal to zero and those for the Kevlar 49-epoxy 
material are negative and have much larger magnitudes 
than the corresponding values for the other materials. 

The results in Figs. 26 and 28 indicate that the val- 

ues of ( 1+Cl ^ 1+Tl ) are positive for the [(±45/0 2 )J s lam- 
inates, for all material systems, except the P-100/3502 
pitch-epoxy and Kevlar 49-epoxy materials, which 

means that fti/(l + Ci)> 0 as long as (1 + T,) > 0. Sim- 
ilarly, the results in Figs. 27 and 29 indicate that the val- 
ues of f ! + ^ - 1— are also positive for all the materials 
\\ + C 2 ;fti 

except the P-100/3502 pitch-epoxy and Kevlar 49-epoxy 
materials. Thus, for the majority of the materials, the 
laminates are loaded by biaxial compression when uni- 
formly heated and by biaxial tension when uniformly 
cooled, and as a result can buckle only for uniform heat- 
ing as long as (1 + T,) > 0 (case 1 ). For the P-100/3502 
pitch-epoxy and the Kevlar 49-epoxy materials, the lam- 
inates are loaded by axial tension and transverse com- 
pression when uniformly heated and by axial 
compression and transverse tension when uniformly 
cooled. For this situation, the laminates can buckle for 
uniform heating (case 2) or uniform cooling (case 1) as 
long as (1 + Tj) > 0. 

The remaining parameters needed for a buckling 
analysis of the long, axially stiff [(±45/0 2 ) m ] s laminates 
are the the nondimensional flexural orthotropy parame- 
ter P and the flexural anisotropy parameters y and 5. 
The values of P for the [(±45/0 2 )J s laminates made of 
the nine material systems are identical to the correspond- 
ing values for the [(±45/90 2 ) m ] s laminates that are shown 
in Fig. 22. Moreover, the values of y for the [(±45/0 2 ) m ] s 
laminates are identical to the corresponding values of 5 
that are shown in Fig. 24 for the [(±45/90 2 ) m ] s laminates. 
Furthermore, the values of 8 for the [(±45/0 2 ) m ] s lami- 
nates are identical to the corresponding values of y that 
are shown in Fig. 23 for the [(±45/90 2 ) m ] s laminates. 


rameters 


Parameters for [(±9) J s laminates Values of the pa- 
and 


for [(±0) m ] s bal- 


(1+CjKI+T,) ~~ ( 1 + C 2 )( 1 + T | ) 
anced, angle-ply laminates composed of the 9 material 
systems are presented in Figs. 30-35. The results in Figs. 
30 and 31,32 and 33, and 34 and 35 correspond to values 
of the relative-compliance ratios R, = R 2 = 0, 0.2, and 
0.5, respectively, and are independent of the number of 
laminate plies. For all the results in Figs. 30-35, the ra- 
tios of the coefficients of thermal expansion are equal (T, 

= t 2 ). 

The results show an extremely wide variation in 
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(i+cun+T,) and ITTc^TO with material system and 

with fiber angle 0. The largest variations are exhibited 

b y o+cS+t,) > follow by (1 ^(iTf, ) - for a " values 

of the relative-compliance ratios considered. Moreover, 

(1+cixi+T,) and (i + c^ fi -T fl) are generally positive for 
all the material systems except, for the most part, for the 
P-100/3502 pitch-epoxy and Kevlar 49-epoxy materials, 
which are negative for several values of the fiber angle 0. 
Thus, depending on the material system and fiber angle, 
a laminate may be loaded by various combinations of 
tension and compression and may buckle when subjected 
to uniform heating, cooling, or both. Like the results for 
the transversely stiff [(±45/90 2 ) m ] s laminates, the [(±45/ 
0/90)J s quasi-isotropic laminates, and the axially stiff 
[(±45/0 2 )J s laminates, the results in Figs. 30-35 show 
that there is a very wide range of possibilities available 
for tailoring the thermal buckling characteristics of a 
laminate family. 

Values of the nondimensional orthotropy parame- 
ter p and the nondimensional anisotropy parameters y 
and 8 for the [(±0) m ] s laminates made of the same 9 ma- 
terial systems have also been presented graphically in 
Ref. 6. The values of P in Ref. 6 for the [(±0)J S lami- 
nates are independent of the stacking-sequence number 
m. The results in Ref. 6 show a series of curves for P 
that vary dramatically with the fiber angle 0 and the ma- 
terial system. The largest values of, and greatest varia- 
tions in, p are generally exhibited by the laminates made 
of the P-100/3502 pitch-epoxy material. In contrast, the 
smallest values of, and least variations in, p are gener- 
ally exhibited by the laminates made of the boron-alumi- 
num material. The values of y and 5 for the [(±0) m ] s 
laminates are strongly dependent on the stacking-se- 
quence number m. The results in Ref. 6 for the flexural 
anisotropy parameters y and 8 are for m = 1 , which cor- 
responds to the highest degree of flexural anisotropy for 
these laminates. The graphical results in Ref. 6 show a 
large effect of the fiber angle 0 and the material system 
on the degree of flexural anisotropy for the [±01 s lami- 
nates. Like the parameter p, the largest values of, and 
greatest variations in, y and 8 are generally exhibited 
by the laminates made of the P-100/3502 pitch-epoxy 
material, and the smallest values of, and least variations 
in, y and 8 are generally exhibited by the laminates 
made of the boron-aluminum material. A comparison of 
results for [(±0) 6 ] s and [±0] s laminates made of IM7/ 
5260 graphite-bismaleimide material is presented in Ref. 
7 (Fig. 21) that indicates that the flexural anisotropy of 
the [(±0) 6 ] s laminates is negligible compared to that of 
the corresponding [±0] s laminates. A similar trend is ex- 


pected for the corresponding laminates made from most 
of the other eight material systems. 

Buckling Coefficients and Critical Temperature Change 

In the present study, the approach that is used to 
compute critical temperature changes for symmetric 
laminates restrained by the apparatus shown in Fig. 3, in 
a generic setting, follows the scenarios designated herein 
as cases 1 through 4. Essentially, one first selects a can- 
didate laminate and computes the nondimensional pa- 
rameters P, y, and 8 (see Eqs. (l)-(4)). Next, the load 
factor L 2 that defines the thermally induced mechanical 
loads is calculated (e.g., by using Eqs. (47) or (50)) for a 
given restraint-medium configuration. This step in- 
volves the calculation of the stiffness-weighted laminate 

thermal-expansion parameters & \ and & 2 , and the calcu- 
lation of the compliance parameters (1 + C,) and 
(1 + C,)/(1 +C 2 ). The parameter ( 1 + C,)/(1 + C 2 ) repre- 
sents the contribution of the elastic restraining medium 
to the thermally induced transverse load (y-coordinate 
direction shown in Fig. 2b). Equivalently, the load factor 
L 2 can be calculated from the appropriate values of 

<5t|/( 1 +Ci) and a 2 f{ 1 + C 2 ). Then, the critical loading 
parameter is determined from figures such as Figs. 22-26 
in Ref. 7, which give the critical (mechanical) loading 

parameter p tr as a function of p, y, 8, L 2 , and the flexural 
boundary conditions. Finally, the critical temperature 
change is calculated by using, for example, Eqs. (49) or 
(52). 

Examples 

To illustrate the use of the results presented previ- 
ously herein with the generic results presented in Ref. 7, 
first consider a simply supported [(±45/90 2 ) 8 ] s laminate 
made of IM7/5260 graphite-bismaleimide material and 
restrained against thermal expansion or contraction by 
the apparatus depicted in Fig. 3. Moreover, consider a 
restraining medium that is define by the relative-compli- 
ance values Rj = R 2 = 0.5 and ratios of the coefficients 
of thermal expansion that are given by Tj = T 2 > -1 and 
whose specific values will be defined later. For these 
values, the axial and transverse compliance of the re- 
straining medium are half of the respective axial and 
transverse compliances of the plate (twice as stiff in each 
direction). From Figs. 22-24, one can estimate that P * 
1.3 and y=8~0. Effectively, this laminate is flexurally 
orthotropic. Similarly, from a laminate analysis code, or 

from Figs. 20 and 21, one could get <5ti = 2.3 x 10‘ 6 /°F, 

and ft 2 /& i= 0.3. In addition, Fig. 8 gives (1 + Cj)(I + 
T,) = 1.6 and Fig. 9 gives (1 +C,)/(! + C 2 ) - 0.6. 
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Combining these results gives 6 i/(l +Ci) = 1.4(1 + T,) 

X 10' 6 /°F and [ 1 ~ > ' Cl )^ 2 = 0.2. For this laminate, the 
\ l +C2JG1 

only destabilizing loads are biaxial compression loads 
that are obtained for uniform heating (case 1). For this 
case, Lj = 1, L 3 = 0, and L 2 = 0.2 is obtained from Eq. 
(47). Note that L 2 = 0.3 for the same plate fully re- 
strained against thermal expansion or contraction. Thus, 
the compliance of the boundary has substantially re- 
duced the thermally induced transverse load. For L, = 1, 

K x = p cr . The value of K x is obtained from Fig. 23 in Ref. 
7 by interpolating the results in the figure for p = 1 and 
P = 1.5. For P = 1 and L 2 = 0.2, Fig. 23 in Ref. 7 gives 
K x = 3.2. Similarly, for p = 1.5 and L 2 = 0.2, Fig. 23 
gives K x = 3.9. Interpolating these values gives K x = 3.6 

for P = 1.3. Substituting K x = 3.6 for p cr and 
<5ti/( 1 + Ci) = 1.4(1 + Tj) X 10 6 f°F into Eq. (49) gives 

= r~% - x 10 6 l, F. Next, let the plate width be 
given by b = 32 in. such that the plate width-to-thickness 
ratio is b/t = 100. For this plate, 0o~ • This result 

indicates that the critical temperature decreases for T, > 
0, which corresponds to the situation where the restrain- 
ing media is expanding and providing additional com- 
pression loads on the laminate. For T, = 0, the relative 
thermal expansion or contraction of the restraining medi- 
um is negligible. For - 1 < T, < 0, the restraining medium 
is contracting relative to the laminate expansion and the 
critical temperature increases. As T ( approaches -1, the 
expansion of the laminate is matched by expansion of the 
restraining medium and buckling cannot occur. More- 
over, for values of -1 < Tj < 0, the possibility of elastic 
buckling may not exist. The corresponding critical tem- 
perature for the plate fully restrained against thermal ex- 
pansion or contraction (Cj = C 2 = 0) and with a thermally 
neutral restraining medium is approximately 1 12°F. 

Next, consider a simply supported [±45/90 2 ] s flex- 
urally anisotropic laminate made of IM7/5260 graphite- 
bismaleimide material and also restrained against ther- 
mal expansion or contraction by the apparatus depicted 
in Fig. 3 with R, = R 2 = 0.5. Likewise, the ratios of the 
coefficients of thermal expansion that are given by Tj = 
T 2 > - 1 . From Figs. 22-24, one can estimate that p ~ 2.0 
and y = 5 ~ 0.2. From a laminate analysis code, or from 
Figs. 20 and 21, one could get 6c 1 = 2.4 x 10 6 /°F, and 

&2/G1 = 0.5. In addition, Fig. 8 gives (1 + C } )(1 + T,) = 
1 .6 and Fig. 9 gives ( 1 +C l )/(\ + C 2 ) = 0.6. Combin- 
ing these results gives &i/(l + C 1 ) = 1.5(1 + T,) x 1 0 6 / 


°F and f 1 +C - )— == 0.3. For this laminate, the only de- 
V l + C2 

stabilizing loads are again biaxial compression loads 
that are obtained for uniform heating (case 1). For this 
case, L, = 1, L 3 = 0, and L 2 = 0.3 is obtained from Eq. 
(47). Note once again that L 2 = 0.5 for the same plate ful- 
ly restrained against thermal expansion or contraction. 
Thus, the compliance of the boundary has substantially 
reduced the thermally induced transverse load. Again, 

for Lj = 1 , K x = p cr . The value of K x is obtained directly 
from Fig. 30 in Ref. 7 and is given by K x = 3.3, which 
corresponds to a wide-column buckling mode. Substi- 
tuting K x = 3.3 for per and &i/(l + Ci) - 1.5(1 + x 

10' 6 /°F into Eq. (49) gives ^-0S r = x 

Next, let the plate width be given by b = 32 in. such that 
the plate width-to-thickness ratio is b/t = 100. For this 

flexurally anisotropic plate, 0o~ j* . The corre- 
sponding critical temperature for the plate fully re- 
strained against thermal expansion or contraction and 
with a thermally neutral restraining medium is approxi- 
mately 66°F. 

Concluding Remarks 

An analytical approach for synthesizing buckling 
results and behavior for long, balanced and unbalanced 
symmetric laminates that are subjected to uniform heat- 
ing or cooling, and elastically restrained thermal expan- 
sion or contraction has been presented. A 
nondimensional buckling analysis for long flexurally 
anisotropic plates that are subjected to combined loads 
has been described and useful nondimensional parame- 
ters have been presented. In particular, stiffness-weight- 
ed thermal-expansion parameters and nondimensional 
compliance parameters have been presented that can be 
used to determine critical temperatures, for a wide range 
of laminate constructions and restraining media, in terms 
of physically intuitive, well-known mechanical-buckling 
coefficients. Moreover, the effects of membrane orthot- 
ropy and membrane anisotropy on the thermally induced 
prebuckling stress state have been determined for a gen- 
eral elastic restraining medium. 

A large number of results have been presented 
herein for some common laminates that are intended, to 
some extent, to facilitate a structural designer’ s transition 
to the use of the generic buckling design curves that are 
included and discussed in the paper. Many of the results 
were previously unknown. In addition, several results 
have been presented that show the effect of laminate con- 
struction and restraining-media compliance on the buck- 
ling behavior, and several cases are presented that 
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indicate when a laminate will buckle because of uniform 
cooling. Results of this type could be important in the 
design of vehicles that use liquid fuels or cruise at high 
altitudes. The use of generic buckling design curves that 
provide physical insight into the buckling problem of the 
present paper in addition to providing useful design data 
has also been discussed. In addition, examples have 
been presented that demonstrate the use of the generic 
design curves. Overall, the analysis approach and ge- 
neric results that have been presented indentify the ef- 
fects or characteristics of elastically restrained laminate 
thermal expansion or contraction, membrane orthotropy 
and anisotropy, and flexural orthotropy and anisotropy 
on laminated-plate buckling in a very general and unify- 
ing manner. Although the results are based on classical 
laminated-plate theory and have been demonstrated for 
infinitely long plates, the approach is applicable to more 
sophisticated plate theories that incorporate effects such 
as transverse-shear flexibility and can be used for finite- 
length plates in many cases. 
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Table 1. Lamina Properties 


Lamina 

property* 

Material Systems 

Boron 

-A1 

S-glass- 

epoxy 

Kevlar 49- 
epoxy 

IM7/ 

5260 

AS4/ 

3502 

AS4/ 

3501-6 

Boron- 

epoxy 

IM7/ 

PETI-5 

P-100/ 

3502 

El, Msi 

33 

mm 

11.02 




29.58 

20.35 

53.5 

Ep Msi 

21 

1.7 

0.8 

1.457 

1.64 

1.30 

2.68 


0.73 

v LT 

0.23 

0.25 

0.34 

0.258 

0.30 

0.30 

0.23 

0.29 

0.31 

G lt , Msi 

7.0 

0.80 

0.33 

0.860 

0.87 

1.03 

0.81 

0.61 

0.76 

x 10VF 

3.2 

3.5 

-2.22 

0.0125 

0.25 

-0.167 

3.38 

-0.14 

-0.64 

ctj x lO'T'F 

11.0 

11.0 

43.89 

14.91 

16.2 

15.6 

16.83 

16.85 

17.2 


* The symbols L and T denote the longitudinal fiber and transverse matrix directions of a specially orthotropic lamina, respectively. 
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Fig. 1 Sign convention for positive-valued stress resultants. 
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(b) Positive values for equivalent mechanical loads 


Fig. 2 Mechanical loads in a plate elastically restrained against thermal expansion or contraction 
caused by uniform heating or cooling. 
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Fig. 3 Central plate elastically restrained by edge plates against thermal expansion or contraction 
caused by uniform heating or cooling in all plates. 
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Fig. 4 Nondimensional compliance coefficient (1 + Cj)(l + for [(±45/0/90) m ] s laminates made of IM7/5260 
material and elastically restrained from thermal expansion or contraction by edge plates (T, = T 2 , 
m = 1, 2, ... ). 



Fig. 5 Nondimensional compliance coefficient (1 + Cj)(l + T,) for [(±45/90 2 ) m ] s laminates made of IM7/5260 
material and elastically restrained from thermal expansion or contraction by edge plates 
(T, = T 2 , m = 1, 2, ... ). 
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Fig. 6 Nondimensional compliance coefficient (1 + C,)(l + T,) for [(±45/02) m ] s laminates made of IM7/5260 
material and elastically restrained from thermal expansion or contraction by edge plates 

(T, = T 2 , m= 1, 2, ...). 
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Fig. 7 Nondimensional compliance coefficient ratio (1 + C, )/(l + C 2 ) for [(±45/0/90) m ] s , [(±45/0 2 ) m ] I , and 
[(±45/90 2 )J, laminates made of IM7/5260 material and elastically restrained from thermal expansion 
or contraction by by edge plates (T, = T 2 and m = 1, 2, ... ). 


25 

American Institute of Aeronautics and Astronautics 






Nondimensional 

compliance 

coefficient, 

(1 +C,)(1 + T,) 



Fig. 8 Nondimensional compliance coefficient (1 + C,)(l + T,) for [(±45/0/90) m ] s , [(±45/0 2 ) m ] s , and [(±45/90 2 ) m ] s 
laminates made of IM7/5260 material and elastically restrained from thermal expansion 
or contraction by edge plates (R 2 = R,, T, = T 2 , and m = 1, 2, ... ). 
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Fig. 9 Nondimensional compliance coefficient ratio (1 + C, )/(l + C 2 ) for [(±45/0/90) m ] s , [(±45/0,) r J s , and 
[(±45/90 2 ) m ] J laminates made of IM7/5260 material and elastically restrained from thermal expansion 
or contraction by edge plates (R 2 = R,, T, = T 2 , and m = 1, 2, ... ). 
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Fig. 10 Stiffness-weighted laminate thermal-expansion parameter + CiJO + T i)] for [(±©) m ] s 

balanced, angle-ply laminates made of IM7/5260 material and elastically restrained from thermal 
expansion or contraction by edge plates (R 2 = R,, T, = T 2 , m = 1, 2, ...). 
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Fig. 11 Stiffness-weighted laminate thermal-expansion parameter + c 2)( 1 + T i)] for [(±9)J S 
balanced, angle-ply laminates made of IM7/5260 material and elastically restrained from thermal 
expansion or contraction by edge plates (R 2 = R„ T, = T 2 , m = 1, 2, ...)• 
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Fig, 12 Critical temperature change for simply supported and clamped laminates made of IM7/5260 
material, fully restrained against thermal expansion, and subjected to uniform heating. 



Fig. 13 Critical temperature change for simply supported and clamped angle-ply laminates made of IM7/5260 
material, fully restrained against thermal expansion, and subjected to uniform heating (m = 1, 2, ...). 


28 

American Institute of Aeronautics and Astronautics 



Critical- 
temperature- 
change ratio, 

eg r 

©r(i+ c f ) 



Fig. 14 Critical-temperature-change ratio as a function of the nondimensional compliance ratio 
(1 + C, )/(l + C 2 ) for [(±45/0/90) m l J laminates made of IM7/5260 material and elastically restrained 
from thermal expansion or contraction by edge plates. 
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Fig. 15 Critical-temperature-change ratio as a function of the nondimensional compliance ratio 
(1 + C. )/(l + C 2 ) for [(±45/0 2 ) m ] s laminates made of IM7/5260 material and elastically restrained 
from thermal expansion or contraction by edge plates. 
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Fig. 16 Critical-temperature-change ratio as a function of the nondimensional compliance ratio 
(1 + Cj )/(l + C 2 ) for f(±45/90 2 ) m ] s laminates made of IM7/5260 material and elastically restrained 
from thermal expansion or contraction by edge plates. 



Fig. 17 Critical-temperature-change ratio as a function of the nondimensional compliance ratio 
(1 + C, )/(l + C 2 ) for [(±0) m ] s laminates (m > 5) made of IM7/5260 material and elastically restrained 
from thermal expansion or contraction by edge plates. 
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Fig. 18 Effects of lamina material properties on nondimensional compliance coefficient (1 + C,)(l + T,) 
for [( ±45/90 2 ) m ] s laminates elastically restrained from thermal expansion or contraction by edge plates 
(R 2 = R„ T, = T 2 , and m = 1, 2, ... ). 



Fig. 19a Effects of lamina material properties on nondimensional compliance coefficient ratio (1 + C,)/(l + C 2 ) 
for [(±45/90 2 ) m l s laminates elastically restrained from thermal expansion or contraction by edge plates 
(R 2 = R„T 1 = T 2 , and m = 1,2, ... ). 
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Fig. 19b Effects of lamina material properties on nondimensional compliance coefficient ratio (1 + C,)/(l + C 2 ) 
for [(±45/90 2 ) ] laminates elastically restrained from thermal expansion or contraction by edge plates 
(R 2 = R„ T, ="T 2 , and m = 1, 2, ... ). 



Fig. 20 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter a, 
for [(±45/90 2 ) m ] s laminates. 
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Fig. 21a Effects of lamina material properties on parameter ratio 62/61 for [(±45/90 2 ) m ] s laminates. 
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Fig. 21b Effects of lamina material properties on parameter ratio 62/01 for [(±45/90 2 ) m ] s laminates. 
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Fig. 22 Effects of lamina material properties on nondimensiona! flexural orthotropy parameter p for [(±45/90 2 )J s 
and [(±45/0 2 )J s laminates. 



Fig. 23 Effects of lamina material properties on nondimensional flexural anisotropy parameter y for [(±45/90 2 ) J s 
laminates and 5 for f(±45/0 2 ) m ] s laminates. 
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Fig. 24 Effects of lamina material properties on nondimensional flexural anisotropy parameter 5 for [(±45/90 2 )J s 
laminates and y for [(±45/0 2 ) m ] s laminates. 



Fig. 25 Effects of lamina material properties on nondimensional compliance coefficient (1 + C,)(l + T,) 
for [ (±45/0/90) m ] s laminates elastically restrained from thermal expansion or contraction by edge plates 
(R z = R p T, = T 2 , and m = 1, 2, ... ). 
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Fig. 26 Effects of lamina material properties on nondimensional compliance coefficient (1 + C,)(l + T,) 
for [(±45/0j) J laminates elastically restrained from thermal expansion or contraction by edge plates 
(R 2 = R„ T, = and m = 1, 2, ... ). 
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Fig. 27 Effects of lamina material properties on nondimensional compliance coefficient ratio (1 + C,)/(l + C 2 ) 
for [(±45/0 2 )J s laminates elastically restrained from thermal expansion or contraction by edge plates 
(R 2 = R„ T, = T 2 , and m = 1, 2, ... ). 
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Fig. 28 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter 61, 
for [(±45/0 2 ) m ] s laminates. 



Fig. 29a Effects of lamina material properties on parameter ratio 62 /$ i for [(±45/0 2 ) m ] s laminates. 
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Fig. 29b Effects of lamina material properties on parameter ratio ftj/fti for [(±45/0 2 ) m ] s laminates. 
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Fig. 30 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter 
«.|/[(1 + C i)(1 +Ti)] for [(±0) m ] J balanced, angle-ply laminates (R 2 = R, = 0,T, =T r m = 1,2,...). 
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Fig. 31 Effects of lamina material properties on stiffness- weighted laminate thermal-expansion parameter ^ 
a 2 /[(i ♦ C 2 )(1 + Ti)] for [(±0) m ] 5 balanced, angle-ply laminates (R 2 = R, = 0, = T 2 , m =1,2, ...). 



Fiber angle, 6, deg 

Fig. 32 Effects of lamina material properties on stiffness- weighted laminate thermal-expansion parameter 
ai/[(l + C i)(1 + Ti)] f or [(±9)J s balanced, angle-ply laminates (R 2 = R, = 0.2, T, =T 2 , m = 1, 2, ...). 
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Fig. 33 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter 
a 2 /[(l + C 2 )(1 + Ti)] for [(±6) m ] 5 balanced, angle-ply laminates (R 2 = R, = 0.2, T, = T 2 , m =1,2, ...). 
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Fig. 34 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter 
ai/[(l + Ci)(l + Ti )] for [(±0) m ], balanced, angle-ply laminates (R 2 = R, = 0.5, T, = T 2 , m = 1, 2, ...). 
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Fig. 35 Effects of lamina material properties on stiffness-weighted laminate thermal-expansion parameter 
a 2 /[(i + C 2 )(l + Ti)] for [(±e)J, balanced, angle-ply laminates (R 2 = R, = 0.5, T, = T 2 , m = 1, 2, ...). 
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